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Abstract 


A  general  method  for  designing  a  laser  system  is  presented.  Using  the  Ti:sapphire  laser 
as  an  example,  the  requirements  of  stability,  astigmatic  compensation,  and  matching 
of  the  pump  and  cavity  modes  are  addressed.  Investigations  into  the  relaxation 
oscillations  of  a  Ti:sapphire  laser  are  reported.  Using  four  level  laser  rate  equation 
theory,  a  technique  is  developed  for  analyzing  relaxation  oscillations  exhibited  by  a 
laser.  This  technique  presents  a  new  and  simple  method  for  measuring  the  upper 
state  lifetime  and  intrinsic  cavity  losses  of  a  laser  system.  Beam-like  vector  solutions 
to  Maxwell’s  equations  are  also  presented.  These  solutions  present  a  more  detailed 
description  of  the  polarization  properties  of  laser  beams.  Experimental  evidence  of 
these  properties  is  shown  using  an  Argon  laser. 


SYSTEM  DESIGN  AND  RELAXATION  OSCILLATIONS 
OF  A  TITANH  M-SAPPHIRE  LASER 


SYSTEM  DESIGN  AND  RELAXATION  OSCTLLATIONS 
OF  A  nTANH'M-SAPPHIRE  LASER 


A  thesis  submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 
Master  of  Science 


by 


WILLIAM  L.  ERIKSON.  B.S. 
United  States  Air  Force  Academy,  1982 


August,  1992 
University  of  Arkansas 


This  thesis  is  approved 
for  recommendation  to  the 
Graduate  Council. 


Thesis  Director: 


Dr.  Surendra  P.  Singh 


Thesis  Committee: 


cl 


Dr.  Rajendra  Gupta 


Dr.  Bill  Durham,  Chemistry 


THESIS  OrPLlCATION  RELEASE 


I  hereby  authorize  the  Eniversity  of  Arkansas  Libraries  to  duplicate  this  thesis  vvlien 
needed  for  research  and/or  scholarship. 

.Agre«l _ 

Refused _ 


Acknowledgements 


I  thank  my  wife  Melissa  lor  lu'r  loving  support  and  faithful  enroiiragement  during 
my  graduate  sc'hool  experience.  I  also  owe  many  thanks  to  my  parents,  who  have 
confidently  supported  me  in  every  <hallenge  I  have  und('rtak('n. 

I  am  indehted  to  Dr.  Singh  for  his  leadership  in  my  research  projei  t.  His  example 
of  pursuing  a  deeper  understanding  of  tlie  world  of  lasers  has  inspired  me  to  do  the 
same.  I  would  also  like  to  thank  my  colleague  ^'ujaing  Qu  for  his  friendship  and  eager 
assistance  in  data  collecting  and  solving  the  endless  little  problems  that  d<'V('lop  in 
experimental  work. 

I  his  thesis  was  typeset  by  the  author  using  IAT£X. 


Contents 


Acknowledgements 


IV 


1  Introduction  1 

2  Propagation  of  Laser  Beams  4 

2.1  Propagation  of  Electromagnetic  Waves .  1 

2.2  The  Para.xial  Wave  Equation .  7 

2.3  Gaussian  Beam  Solutions .  !) 

2.4  Transformation  Characteristics  of  Gaussian  Beams  .  13 


3  Laser  Design  16 

3.1  Ring  versus  Standing  Wave  Cavity  .  l(i 

3.2  Stable  Cavity:  The  Ray  Matri.x  .Approach  .  17 

3.3  Intracavity  Beam  Waists:  .Astigmatic  Compensation .  25 

3.4  .Mode  Matching .  30 

3.5  Intracavity  Elements .  .35 

3.5.1  Birefringence  .  35 

3.5.2  1-  nidirectional  Lasing .  30 

3.5.3  Tuning  to  Single  Wavelength .  3S 

4  Laser  Dynamics  41 

1.1  f.aser  Rate  Equations .  11 

1.2  Rela.xation  Oscillations:  Laser  Rate  Equation  .Analysis  .  11 

1.3  Rela.xation  Oscillations:  .Numerical  Modeling .  17 

1.1  .A  Method  for  Determining  Intrinsic  Cavity  Loss .  50 


5  Experimental  Results  53 

').!  Laser  Deseript ion .  l;! 

I'.xpei  iiiient  1:  H('laxation  Oscillatioi.s  .  oo 

Experimental  Eroeednre .  '>() 

Hesiilts .  ')^ 

0.2. 3  Discussion  of  Results .  fio 

0.2.1  Results  of  .Xumeriral  .Moileling  .  LI) 

0.3  Experiment  2;  Cavity  Los.ses .  Tfi 

0.3.1  Experimental  Setup:  .\0.\1  Calibration .  Tti 

0.3.2  Experimental  Setup:  (.’avity  Losses .  77 

0.3.3  Results .  71) 

•0.3.4  Discussion  of  Results .  S2 

6  Polarization  Properties  of  Gaussian  Beams  84 

6.1  Paraxial  Solutions  to  Maxwell's  Equations .  So 

6.2  Linearly  Polarized  Electrontagnetic  Waves .  S7 

6.3  Circularly  Polarized  Electromagnetic  Waves  .  !)2 

6. 1  Experimental  Observations .  t)2 

7  Conclusions  96 

Bibliography  98 


vi 


Chapter  1 


Introduction 


The  Tiisapphire  laser  has  been  the  subject  of  continuing  research  since  it  was  first 
introduced  by  Moulton  [1]  as  a  broadly  tunable  solid  state  laser.  Because  of  it' 
broad  tunability  (ss  650-1200  nm).  high  outptit  power,  and  demonstrated  stabilifc. 
the  Titsapphire  laser  is  gaining  prominence  in  many  areas  of  research.  I  hese  areas 
include  high  resolution  spectroscopy  and  high  energy  laser  applications. 

Many  areas  of  high  resolution  spectroscopy  research  require  frequency  stabilizcnl. 
single-mode,  continuous  wave  (cw')  lasers.  Although  dye  lasers  have  been  dominant 
in  this  field,  the  Tirsap  phi  re  laser  offers  the  advantages  of  wide  tunability.  ease  of  use. 
and  increased  stability.  The  dominant  absorption  band  for  Titsapphire  is  centered 
around  500  nm.  allowing  direct  pumping  with  an  argon  ion  laser  for  cw  operation  and 
doubled  Nd:YAG  lasers,  dye  lasers,  and  fla.sh  lamps  for  pulsed  operation,  .\ctive  fre 
quency  stabilization  involves  a  reference  cavity  (Fabry  Perot),  some  type  of  feedback 
to  the  laser  cavity,  and  an  intracavity  compensating  element.  By  actively  stabilizing 
a  Ti:sapphire  laser,  a  frequency  stability  of  1  kHz  rms  relative  to  the  reference  cavity 
has  been  observed  [2].  Stabilized  Ti:sapphire  lasers  have  been  used  for  examining 
the  hyperfine  structure  of  the  D  lines  and  the  absorption  spectroscopy  of  the  5S— >51) 
transition  in  rubidium  [3|.  .A  similar  Tirsapphire  laser  was  used  to  investigate  the 
two-photon  excitation  of  the  2S— +4S  transition  in  atomic  hydrogen  [4|.  Other  areas 
of  spectroscopic  research  require  broadly  tunable  short  laser  pulses  of  relatively  high 
energy.  This  is  another  area  where  the  capabilities  of  the  Ti:sapphire  laser  are  being 
exploited.  Various  methods  for  active  and  passive  mode  locking  of  a  Ti:sapphire  laser 
have  been  employed  to  generate  high  power  pulses  in  the  femtosecond  regime  [5,  6.  7k 

High  energy  laser  applications,  such  as  future  remote  space  sensing  systems,  have 
generated  the  need  for  studying  the  energy  output  capabilities  of  Titsapphire  la.'iers  [Si. 
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Output  powers  of  0.3  terawatt  for  123- fs  pulses  have  been  observed  using  the  technique 
of  chir[)ed-pulse  amplification  [9],  For  cw  operation.  330  watts  of  output  power  has 
been  proposed  using  a  liquid  nitrogen-cooled  Tirsapphire  laser  [10]. 

Other  areas  of  research  where  Tirsapphire  lasers  may  be  useful  are  nonlinear  of)t  if  s. 
laser  cooling,  and  trapping  of  atoms.  Because  of  the  many  possible  and  diverse' 
applications  for  a  Tirsapphire  laser,  there  is  a  need  for  a  general  method  for  designing 
simple  and  flexible  Tirsapphire  laser  cavities  from  which  a  specific  laser  system  can 
be  constructed. 

Although  most  of  the  current  research  involving  the  Tirsapphire  laser  is  centered 
around  areas  of  specific  application,  little  attention  has  been  paid  to  the  study  of  the 
dynamical  behavior  of  this  laser  system.  One  particular  dynamical  behavior  of  interest 
is  the  relaxation  oscillations  exhibited  by  the  Tirsapphire  laser.  These  oscillations 
characterize  the  behavior  of  certain  laser  systems  when  they  are  perturbed  from 
steady  state.  Relaxation  oscillations  have  been  observed  and  studied  in  a  number  of 
other  laser  systems  including  ruby  [11],  dye  [12]  and  copper- vapor  [13].  By  studying 
these  relaxation  oscillations,  one  can  gain  a  greater  insight  into  the  parameters  that 
drive  the  dynamics  of  the  laser.  Recent  studies  in  the  ultra-high  speed  relaxation 
oscillations  (up  to  39  GHz)  of  vertical  cavity  surface  emitting  diode  lasers  indicate 
potential  application  in  the  area  of  optical  interconnect  [14]. 

In  this  thesis,  w?  present  a  general  method  for  designing  a  Tirsapphire  laser  that 
can  be  customized  for  specific  research  applications.  We  also  report  on  investigations 
of  the  relaxation  oscillations  of  a  Tirsapphire  laser  that  we  constructed.  To  our  knowl¬ 
edge.  these  are  the  first  such  investigations  to  have  been  carried  out  on  a  Tirsapphire 
la.ser  system.  By  analyzing  the  relaxation  oscillations,  we  measured  the  upper  state 
lifetime  and  the  intrinsic  cavity  losses  of  our  laser.  Our  method  for  analyzing  these 
oscillations  presents  a  simple  and  new^  technique  for  measuring  these  important  laser 
parameters. 

In  Chapter  2.  we  review  the  propagation  characteristics  of  laser  beams  including 
their  transformation  characteristics  through  optical  elements.  With  this  foundation, 
we  describe  a  technique  using  ray  transformation  matrices  in  Chapter  3  to  address  the 
problems  of  stability,  astigmatism,  and  mode  matching  in  the  design  of  a  laser  system. 


W’e  also  describe  the  design  of  a  specific  Ti:sapphire  laser.  In  Chapter  4  we  discuss 
the  theory  of  relaxation  oscillations  based  on  two  simplified  laser  rate  etpiations.  W’e 
then  develop  a  method  for  analyzing  the  relaxation  oscillations  exhibited  by  certain 
lasers  which  allows  us  to  measure  the  upper  state  lifetime  and  the  intrinsic  cavity 
losses  of  a  laser  system.  Chapter  o  describes  two  sets  of  experiments  conducted  with 
our  Titsapphire  laser.  The  first  set  involved  the  study  of  the  relaxation  oscillations 
generated  by  cavity  loss  modulation.  The  second  set  of  experiments  was  conducted  to 
determine  the  intrinsic  cavity  losses  of  our  laser  by  an  independent  method.  Finally, 
in  Chapter  6  we  describe  vector  beam-like  solutions  of  Ma.xwell's  equations.  W'e  use 
these  solutions  to  discuss  the  polarization  properties  of  (4aussian  laser  beams.  We 
also  present  experimental  evidence  representing  these  solutions. 
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Chapter  2 


Propagation  of  Laser  Beams 


Any  study  of  lasers  must  involve  a  basic  understanding  of  the  propagation  character¬ 
istics  of  laser  beams.  To  this  end.  we  begin  this  chapter  with  a  review  of  the  wave 
nature  of  light  and  its  propagation  characteristics  in  free  space  (vacuum).  VVe  then 
derive  the  paraxial  wave  equation,  which  describes  propagation  of  electromagnetic 
waves  (beams)  whose  energy  is  concentrated  near  the  axis  of  propagation,  and  look 
at  a  particular  beam-like  solution  to  this  equation.  Finally,  we  will  see  how  laser 
beam  parameters  are  transformed  as  they  pass  through  various  optical  elements. 

2.1  Propagation  of  Electromagnetic  Waves 

The  behavior  of  an  electromagnetic  field  is  governed  by  .Maxwell's  equations,  which 
for  free  space  can  be  written  in  the  form 


VE(r.O  = 

0, 

(•2.1) 

VxE(r,<)  = 

^B(r.  0 
df 

(2.2) 

VB(r,<)  = 

0  . 

(•2.;4) 

VxB(r./)  = 

1  dEir.f) 
dt 

(2.4) 

where  c  is  the  speed  of  light  in  free  space.  E(r.  #)  is  the  electric  field  and  B(r.  0  is 
the  magnetic  field.  If  we  eliminate  B(r,  0  from  these  equations  by  taking  the  curl 
of  Fq.  (2.2)  and  the  time  derivative  of  Eq.  (2.4).  we  obtain  a  closed  equation  for  Uie 
electric  field.  Using  the  fact  that  time  and  space  derivatives  commute  and  the  vector 
identity 

V  X  (V  X  E)  =  V(V-E)  -  V^E.  (2..-)) 
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togclhor  with  Fwi.  (2.1),  vve  find  that  the  electric  field  satisfies  the  irai'f  fqiiaitoit 

Similarly,  by  taking  the  ctirl  of  Eq.  (2.4)  and  the  time  derivative  of  E([.  (2.2)  we 
can  show  that  the  magnetic  field  also  satisfies  the  wave  ecjnation.  This  means  that 
Maxwell's  ecjuations  admit  wave-like  solutions  and  each  component  of  the  E  field  and 
B  field  satisfies  the  scalar  wave  equation 

=  0.  (2.71 

The  solutions  of  this  equation  are  of  the  form: 

{f  {t  z/c)  plane  wave 

spherical  wave  (2.8) 

T  pje)  cylindrical  wave 

where  'f(r,  t)  represents  any  component  of  E(r.  tl)  or  B(r, /)•  The  particular  form  of 
the  solution  chosen  depends  on  the  symmetry  of  the  problem.  In  describing  laser 
propagation,  we  are  interested  in  quasi-monochromatic  fields  (coherent  light),  so  we 
desire  solutions  of  the  form 

'l'(r,0  =  (2.!J) 

If  we  substitute  this  solution  into  the  scalar  wave  equation  (2.7),  we  find  the  space 
dependent  part.  'Pc(r),  satisfies  the  Helmholtz  equation 

(v2  +  :^)vlr„(r)  =  0.  (2.10) 

Perhaps  the  most  familiar  solution  of  this  equation  is  the  plane  wave  solution.  VVe 
can  represent  a  plane  wave  propagating  in  the  r-direction  by 

'{'.(r)  =  eve’'''' ,  (2.11) 

where  v.j  is  a  constant  and  the  propagation  constant  k  is  related  to  the  wavelength 
and  angular  frequency  of  the  wave  by 
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The  E  field  of  a  plane  electroniagnetie  wave  is  then  written  as 


E(r.O  =  . 


(2. Id) 


and  the  B  field  as 

B(r./)  =  A-  X  0  .  (2,1  n 

c 

For  plane  waves.  Eq  is  a  constant  vector.  These  two  equations  satisfy  F.qs.  (2.2)  and 
(2.4).  Ma.xwell’s  equations  (2.1)  and  (2.3t  require 


i-  ■  E(r,0  =  A-  •  B(r. 0  =  0  . 


(2.15) 


Then  the  energy  flow  as  defined  by  the  Poynting  vector 

S  =  -^?Re(E*  X  B)  =  \eocEfk  (2.16) 

2fio  2 

is  in  the  ^-direction  as  expected,  .■\lthough  the  unidirectional  propagation  of  plane 
waves  is  characteristic  of  laser  beams,  the  fact  that  the  wave  has  the  same  amplitude. 
'l'o(r)  =  for  the  entire  plane  s  =  constant  is  not  characteristic  of  a  laser  beam. 

In  addition  to  a  predominant  direction  of  propagation,  laser  beams  also  have  a  finite 
extent  in  the  transverse  direction.  VVe  therefore  look  at  other  potential  solutions  of 
Eq.  (2.10)  that  mimic  these  properties  of  laser  beams. 

The  second  important  solution  of  Eq.  (2.10)  is  the  spherical  wave,  which  can  be 
written  for  r  ^  0  in  the  form 

»I>,(r)  = (2.17) 
r 

where  /I  is  a  constant.  This  represents  a  wave  of  constant  amplitude  over  a  sphere 
of  radius  r.  Since  we  are  interested  in  the  predominantly  unidirectional  propagation 
of  a  wave,  we  restrict  ourselves  to  a  small  cross-section  of  the  wave  near  the  r-axis 
(paraxial)  at  a  distance  z  =  R  away  from  the  origin,  or  the  point  source.  For  these 
points  close  to  the  r-axis,  we  can  write  r  as 


r  =  R 


£i±J!!V" 

) 


(2.18) 


and  since  -f  R^  for  the  para.xial  points,  we  can  approximate  r.  using  the 

binomial  expansion,  by 


r  %  /?  + 


2R 


(2.19) 
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I  he  N()heriral  wave  solution  on  the  c-axis  at  :  =  R  can  then  be  written  as 

=  ■^f<kR^tk(T‘+y^)l2R  (2. 20) 

R 

for  distances  far  from  the  source.  Because  of  its  small  relative  magnitude  compared  to 
the  second  term  in  Eq.  (2.19)  can  be  neglected  in  the  denominator  of  Eq.  (2.20). 
It  must  be  kept  in  the  exponential  term  however,  because  there  it  is  compared  to 
a  wavelength.  In  Eq.  (2.20),  R  is  the  radius  of  curvature  of  the  phase  fronts.  We 
now  have  a  solution  to  the  Helmholtz  equation  which  has  a  non-uniform  amplitude 
in  a  plane  perpendicular  to  the  r-axis.  We  will  see  that  this  solution  suggests  very 
important  "beam-like"  solutions  to  the  paraxial  wave  equation,  which  we  derive  next. 


2.2  The  Paraxial  Wave  Equation 

We  are  interested  in  solutions  that  give  a  finite  transverse  extent  and  yet  travel 
predominantly  in  the  ^-direction.  To  this  end  we  propose  a  solution  of  the  form 

<p,(r)  =  t>(r)e‘^''  (2.21) 


where  c’(r)  describes  the  transverse  profile  of  the  beam.  Substituting  this  into  the 
Helmholtz  equation  (2.10)  and  using  the  fact  that  k  =  ^'/c.  we  obtain 


t!'(r)  +  2ik 


dti'(r) 

dz 


+ 


d^0{T) 

dz^ 


ikz 


=  0 . 


(2.22) 


We  now  assume  that  the  c  dependence  of  ii'{r)  is  slow.  This  means  that  the  transverse 
profile  of  the  beam  does  not  vary  significantly  over  distances  comparable  to  the  optical 
wavelength  A  =  'Iir/k.  In  other  words,  the  beam  spreads  slowly  as  it  propagates  in 
the  c  direction.  This  assumption  is  known  as  the  paraxial  approximation  and  is 
represented  mathematically  by  the  conditions 


(?u'(r) 

dz 

d^n'(r) 

dz^ 


<  *:|c(r)|  . 
I^v(r) 


<  2k 


and 


Vff(r) 


(2.2;i) 

(2.24) 
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wlit’n*  tti«'  traiisrtrsf  Laplacian  is  given  In 


if  if 
i).r^  ^  i)!f  ■ 


( '2.2') ) 


rnder  t  lie  paraxial  approximat  ion.  t  lie  llelmlioll  z  «’(|iiat  ion  (2.10)  leads  to  t  lu'  ixiid.nitl 
wart  f(f nation: 

,  i)i.'(T) 

VTi’(r)  +  2ik — ; —  =  0  .  (2.2()) 

i)z 

Before  wo  look  at  specific  solutions  to  this  e<iuation.  let  us  invi'st igate  tin'  validity 
of  the  paraxial  wave  approximation.  Sipipose  we  have  a  wave  traviding  in  the 
plane  whose  A'-vector  makes  an  angle  B  with  the  r-axis.  Since  we  can  always  descrilie 
a  wave  as  a  superposition  of  plane  waves,  we  can  represent  the  wave  component  in 
the  x-z  plane  as 

where 


VV’e  then  find  for  each  term  in  Eq.  (2.24): 


ifv{x.  c) 
dz^ 

Uz 

z) 

ilr^ 


-k-^ico.sO  -  D^dJ-.r) 
— 2A'^(co,s^  —  1  ') 


(2.20) 

(2.40) 

(2.41) 


I  sing  the  approximations  sinO  ss  0  and  co.'^O  —  1  ^  0^/2.  we  find  that  Lfp  (2.21)  leads 
to  t he  ineipiality 

I  hus  the  paraxial  apjiroximation  is  good  as  long  as  this  ine(|uality  is  satisfied.  It 
we  say  that  two  orders  of  magnitude  is  sufficient  to  call  the  approximation  valid, 
then  the  ineijuality  holds  for  0  <  0.2  railian  or  12“^.  This  means  the  paraxial  wavi' 
approximation  is  valid  as  long  as  the  beam  does  not  converge  or  diverge  outside  of  a 
cone  of  %  2T'.  or  as  long  as  most  of  the  plane  waves  comprising  the  beam  hav('  their 
A-V(x'tor  inside  a  cone  of  %  2f". 
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2.3  Gaussian  Beam  Solutions 


\\V  now  recall  the  spherical  wave  solution.  Kq.  (2.20)  derived  in  Section  2.1.  Since  it 
is  a  solution  to  the  Helmholtz  equation  (2.10)  in  the  para.xial  appro.\iination.  it  must 
also  satisfy  the  paraxial  wave  equation.  VVe  generalize  the  form  of  the  spherical  wave 
solution  such  that  it  has  a  Gaussian  transverse  profile. 

^  (2.33) 


where  A  is  a  constant  and  we  introduce  the  complex  beam  parameter  f/(~)  and  a 
complex  phase  shift  p{z).  Requiring  Eq.  (2.33)  to  be  a  solution  of  Eq.  (2.26)  alluw.s 
us  to  determine  the  parameters  q(:)  and  p(c).  With  i!’(r)  given  by  Eq.  (2. .33).  each 
term  in  Eq.  (2.26)  can  be  written  as 


dx^ 


«/’(r)  = 


ik 


k^x^ 


|v(r)  = 


.<?(>) 

'  ik  k^y^  ^ 


t-’(r) 

e(r) 


,,  (j2  -p  y^)dq{z)  ,  .dp{z)\ 
-tk—r-rr- - : —  +  i— —  1  t'lr) 


2qHz}  dz 


dz 


(•2.34) 

(2..3.5) 

(2.36) 


With  the  help  of  these,  we  find  that  Eq.  (2.26)  becomes 


'(iqjz) 

dz 


-  I  -  2it 


'dpjz) 

dz 


To  satisfy  this  equation,  we  must  simultaneously  have 


dq{:] 

dz 


=  1 


and 


dp(  =  ) 
dz 


q(=) 


Integrating  Eq.  (2.38),  we  find 


?(-)  =  qo  +  = 
qo  +  - 


p(z)  =  I  In 


q^ 


where  q{0)  =  q^ 
where  p(0)  =  0  . 


(2.37) 


(2.38) 


(2.39) 

(2.40) 


In  order  to  understand  the  physical  significance  of  the  complex  beam  parameter 
q{z).  vve  write  it  in  terms  of  its  real  and  imaginary  parts  by  introducing  two  real  beam 


9 


paraiiiPtf^r'i  Hi:)  and  »'(-) 


1  _  I  .  A 

<7(r)  R{:)  ^  '  KU’^:)  ' 

In  the  plane  c  =  0  this  equation  leads  to 

<7(0 )  Ro'^^TTwy 


(2.  U  ) 


(1^-12) 


If  we  choose  r  =  0  to  be  the  point  where  Ro  =  3C.  then  q,,  is  pure  imaginary  and  we 
can  write 


where 


1  - 

<io  Jric2 

-  _  £!ii  _ 

A  ■ 


•>  ■ 


(2.  Id) 


(2.14) 


Using  Eqs.  (2.43)  and  (2.44)  we  can  rewrite  Eq.  (2.39)  &s 


q{z)  =  -izo  +  c 


(2,4o) 


and  Eq.  (2.41)  as 


1 


1 


l/9o 


qi')  9o  +  c  l  +  c/g^' 

Separating  the  real  and  imaginary  parts  of  Eq.  (2.46)  we  find 


1 


+  i- 


1/3. 


<7(-)  1+cV.i  1  +  --V 

Comparing  Eq.  (2.47)  with  Eq.  (2.41)  we  find  that 

,2 


2/^2 

o 


/?(c)  =  c  + 


and 


u’U  =  «• 


1 


-3). 


(•2.I6) 


(2.47) 


(2.48) 


(2.19) 


Making  use  of  Erp  (2.43)  in  Eq.  (2.40)  vve  can  write  the  expression  for  p{z)  as 
p(  c  )  =  /  In  =  i  In 


I  +  —  -  tan' 


lO 


(2.A0) 


Finally,  using  Fqs.  (2.41)  and  (2.50)  in  Eq.  (2.3:4)  and  letting  4-  .t/"’.  we  ran 

write  tlie  (laussian  beam  solution  as 


r  r  = 


+  -'V-i 


(2.51) 


In  terms  of  the  expression  for  u'{:)  from  Eq.  (2.49)  we  write  the  full  solution  to  the 
paraxial  wave  equation  in  the  form 

^^(r)  =  *■!'?  g  pV  ^  tan  r., )]  ^ ^ 

(/•(;) 


We  are  now  in  a  position  to  understand  the  physical  signifirance  of  w(:).  .  and 

R{z)  and  how  they  characterize  Gaussian  beams.  If  we  look  at  the  intensity  of  a  wave 
of  total  power  f*, 


;(r)  =  -«„|£(r)l"  = 


~2p^/w^{z) 


wHz) 


T^IL-^Z) 


(2.53) 


we  see  a  Gaussian  distribution  in  the  transverse  direction  with  jrf:-)  representing  the 
characteristic  width  of  the  beam  or  the  beam's  spot  size  at  the  plane  intersecting  the 
beam  a.xis  at  c.  The  meaning  of  "spot  size"  is  not  uniform  in  the  literature.  F'or 
our  purposes,  we  define  the  spot  size  to  be  ir(c)  and  therefore  2w{z)  represents  the 
diameter  of  the  beam.  The  plane  intersecting  the  beam  axis  at  r  =  0  is  known  as  the 
beam  waist.  From  Eq.  (2.49)  we  see  that  the  spot  size  is  minimum  at  the  beam  waist 


w{Q)  =  Wo  .  ('2.54) 

.\  diagram  of  the  beam  waist  region  of  a  Gaussian  beam  is  shown  in  Fig.  2.1.  The 
length  scale  over  which  the  beam's  spot  size  grows  in  the  direction  of  propagation  is 
determined  by  the  quantity  Zg,  known  as  the  Rayleigh  range,  defined  by  Eq.  (2.44). 
From  Eq.  (2.49)  we  can  see  that 

it'izo}  =  \/2wo  (2.55) 

so  the  cross  section  of  the  beam  doubles  over  the  length  equal  to  the  Rayleigh  range. 
The  distance  of  2zo  from  — to  Zg  is  referred  to  as  the  confocal  parameter  of  the 
beam,  the  distance  over  which  the  beam  remains  approximately  collimated. 
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Figure  2.1:  A  diagram  of  the  beam  waist  region  of  a  Gaussian  beam.  The  beam  waist 
is  the  plane  r  =  0.  The  minimum  spot  size  is  found  at  the  beam  waist  and  has  a 
diameter  of  2u’<,.  The  angle  of  divergence  of  the  beam  is  represented  by  0  and  the 
radius  of  curvature  of  the  phase  fronts  is  given  by  R(z). 


To  understand  the  physical  significance  of  R(z).  we  compare  the  full  Gaussian 
beam  solution.  Eq.  (2.52).  with  the  spherical  wave  solution.  Eq.  (2.20)  and  see  that 
R{z)  represents  the  radius  of  curvature  of  the  phase  front.  This  radius  of  curvature  i.s 
infinite  at  the  beam  waist,  where  the  beam  is  collimated.  In  the  far  field,  the  radius 
of  curvature  is  approximately  equal  to  c.  By  convention,  for  a  beam  traveling  in  the 
positive  c  direction.  R{  z)  >  0  for  a  diverging  beam  and  R(z)  <  0  for  a  converging 
beam.  It  is  interesting  to  note  that  if  the  quantity  and  the  beam  waist  location 
is  known,  all  other  information  concerning  the  Gaussian  beam  at  any  other  point  in 
space  can  be  determined  from  Eqs.  (2.44)  -  (2.49).  This  means  Zo  uniquely  determines 
a  Gaussian  beam. 

We  saw  in  Section  2.2  that  the  paraxial  approximation  would  be  valid  if  the 
solution  did  not  diverge  at  an  angle  greater  than  approximately  Ti'"  (or  %  1  /4  radian). 
VVe  solve  for  the  beam  divergence  angle  of  the  Gaussian  spherical  wave  solution  by 
evaluating  Eq.  (2.49)  at  c  ^  z^.  With  the  help  of  the  binomial  expansion  we  find 


u-(-'  >  Co)  = 

The  angle  of  divergence  is  expressed  as 

0  ss  tan  9  = 


aiz) 


zA 

A 


(2.56) 


(2.57) 
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I  bus  as  long  as  the  minimum  spot  size  is  larger  than  a  few  wavelengths,  a  condition 
which  is  fulfilled  for  most  laser  beams,  the  paraxial  approximation  is  expected  to 
hold. 

1  he  solution  we  have  assumed  in  this  section  is  the  so  called  fundamental  solu¬ 
tion.  In  general,  there  are  other  solutions  of  the  paraxial  wave  equation  (2.J6).  For 
rectangular  s\’nunetry.  there  are  the  Hermite-Gaussian  solutions,  which  have  the  form 

'-'mn(r)  =  f  f  exp[-dm  -h  U  +  1 )/«  U  “ '  (  C  /  C  j]  X 

exp[ik{x^  -f  i/)/-2R{z)  -  (j-^  +  (/^)/ic^(c)]  .  (2..'),'^) 

Here  Hni-r)  represents  a  Hermite  polynomial  of  order  n  and  argument  x.  For  cylin¬ 
drical  symmetry,  the  solutions  are  known  as  Laguerre-Gaussian  solutions  and  are  of 
the  form 


X  e\p[ikp^ /2R(z)  —  p^fiu^{z)  +  imo]  .  (2.')!f) 

where  /./"*(  j)  refers  to  an  associated  Laguerre  polynomial  of  degree  /.  For  most  appli¬ 
cations,  the  fundamental  solution  is  the  desired  solution,  and  is  therefore  most  often 
encountered  in  the  literature.  We  now  look  at  the  effect  of  various  optical  elements 
on  a  Gaussian  beam. 

2.4  Transformation  Characteristics  of  Gaussian  Beams 

V\e  have  seen  that  for  free  space  propagation  laser  beams  may  be  represented  by  a 
Gaussian  spherical  beam  solution  to  the  paraxial  wave  equation.  We  now  investigate 
how  these  beams  are  transformed  as  they  pass  through  various  optical  elements. 
.Specifically,  we  are  interested  in  how  q{z)  transforms  as  the  beam  propagates. 

W'e  start  with  the  example  of  free  space  propagation.  .According  to  Eq.  (2.39). 
an  initial  beam  parameter,  q,.  will  transform  into  the  final  beam  parameter  qj  after 
propagating  through  a  distance  d  by 

<}f  =  A  d  .  (2.fi0) 
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If  UT  represent  this  transformation  by  the  equation 


'// 


-H  B 
f  (h  +  ^ 


we  can  write  the  coefficients  in  the  beam  transformation  matrix 


(2.bl  ) 


’  .4 

B  ' 

'  \  cl' 

C 

D. 

0  1 

We  now  take  the  example  of  a  thin  lens  with  a  focal  length  /.  which  transforms  a 
spherical  wave  of  radius  /?,  incident  from  the  left  into  a  spherical  wave  of  radius  Rj 
exiting  the  lens  to  the  right  according  to  the  equation 


_L  _  _L  _  1 

Rj~  R,~  /■ 


For  a  Gaussian  beam  incident  on  a  thin  lens,  we  expect  the  radius  of  curvature  to 
be  transformed  as  in  Eq.  (2.63).  but  the  spot  size  just  to  the  left  of  the  lens  should 
be  exactly  the  same  as  the  spot  size  just  to  the  right  of  the  lens.  We  can  therefore 
relate  the  transformed  beam  parameter,  9/,  just  after  the  lens  to  the  incident  beam 
parameter,  q-,.  just  before  the  lens  by 


_L  -  1  _  1 

9/  9.  /  ' 

Solving  explicitly  for  q/, 

= 

f  ~  Qi  -Qt/ f  +  ^ 

we  can  write  the  beam  transformation  matrix  for  a  thin  lens 


'  A  b' 

1  0 ' 

r  D 

.  '1//  ‘  . 

(2.64) 


(2.63) 


(2.66) 


These  two  specific  transformations  are  examples  of  the  general  ABC D  law  for 
Gaussian  beams  [1.5].  The  transformed  q  parameter  of  a  Gaussian  beam  can  be 
obtained  from  Eq.  (2.61)  where  the  coefficients  are  determined  from  the  ray  transfer 
matrices  of  geometric  optics. 

In  this  chapter  we  have  reviewed  the  propagation  of  electromagnetic  waves  arul 
saw  how  they  can  be  used  to  represent  laser  beams.  From  a  scalar  treatment  ot 
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Maxu-ell's  equations,  vve  derived  the  para.xial  wave  equation.  Bearn-like  solutions  to 
the  para.xial  wave  equation  characterize  most  of  the  propagation  characterist ic-;  of 
laser  beams,  for  (laussian  beam  sohitions.  the.se  characteristics  are  determined  from 
the  complex  beam  parameter  q(z).  The  manner  in  which  a  (iaussian  beam  transforms 
through  optical  elements  is  determined  from  the  ray  transfer  matrices  of  gfxmietric 
optics.  These  matrices  will  play  an  important  role  as  we  design  a  laser  system  in  the 
next  chapter. 
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Chapter  3 


Laser  Design 


Solid  state  lasers  have  proved  to  he  widely  tunable  and  simpler  to  use  than  dye 
lasers.  They  are  finding  inrreasing  use  in  spertrosropir  and  nonlinear  optical  experi¬ 
ments.  Because  of  broad  tunability.  in  order  to  minintize  losses,  one  must  use  crystals 
with  Brewster  windows.  This  presents  the  problem  of  how  to  compensate  for  the 
astigmatism  induced  by  Brewster  surfaces  and  other  nonparaxial  elements.  For  most 
applications,  cylindrically  symmetric  (TE.VIqo)  beams  are  needed.  In  this  chapter  we 
present  a  general  method  that  allows  us  to  compensa^''  ""tigmatism  using  curved  mir¬ 
rors.  Our  objective  is  to  design  a  tunable  single -frequency,  single-mode  Titsapphire 
laser.  We  begin  our  discussion  of  laser  system  design  with  a  comparison  belwerui 
standing  and  ring  cavity  designs.  We  then  describe  the  stability  requirements  for 
a  laser  cavity  and  address  astigmatic  compensation.  The  problem  of  matching  the 
pump  mode  into  the  cavity  mode  is  considered  next.  We  use  our  cavity  design  as  an 
example.  Finally,  we  describe  various  intracavity  elements  that  can  be  used  to  make 
the  laser  unidirectional  and  tunable  to  single  wavelength. 

3.1  Ring  versus  Standing  Wave  Cavity 

The  first  consideration  in  designing  a  laser  is  to  choose  between  a  standing  wave  or  a 
ring  cavity  to  optimize  the  particular  laser  system.  .\  standing  wave  cavity  is  generally 
simpler  in  design  and  offers  the  advantage  of  the  electromagnetic  wave  passing  through 
the  gain  medium  twice  during  a  round  trip  in  the  cavity.  One  disadvantage  of  a 
standing  wave  cavity  however  is  spatial  hole  burning.  The  electromagnetic  field  inside 
a  standing  wave  cavity  may  be  considered  to  be  the  superposition  of  two  oppositely 
propagating  traveling  waves.  The  interference  between  these  two  waves  sets  up  nodes 
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inside  the  cavity  (and  hence  inside  the  gain  medium)  where  the  intensity  of  the 
electromagnetic  wave  goes  to  zero.  Between  two  nodes  is  an  antinode  corresponding 
to  intense  fields.  The  result  is  periodic  gain  depletion  or  “holes  burned"  at  the  spacing 
of  A/2  inside  the  gain  medium.  Therefore,  standing  wave  cavities  do  not  utilize  tli(' 
gain  completely  and  various  grating  effects  can  occur.  The  problem  of  spatial  hole 
burning  is  overcome  in  a  ring  laser,  where  no  standing  waves  exist.  Another  advantage 
of  ring  lasers  is  found  in  the  pumping  scheme.  Many  solid  state  lasers  require  optical 
pumping  of  the  gain  medium,  which  is  often  accomplished  with  another  laser.  In  a 
standing  w’ave  cavity,  part  of  the  pump  beam  can  be  reflected  back  onto  itself,  causing 
interference  and  stability  problems,  especially  in  situations  involving  longitudinal 
pumping  schemes.  In  a  ring  laser,  such  reflections  or  feedback  are  not  present.  A 
disadvantage  of  a  ring  laser  is  the  directional  mode  competition  which  may  result 
in  large  intensity  fluctuations.  This  directional  competition  can  be  suppressed  by 
introducing  a  Faraday  rotator  and  waveplate  assembly  as  described  in  Section  3.5.2. 
Finally,  the  design  of  a  ring  l«iser  generally  affords  for  easier  introduction  of  intracavity 
elements  with  the  second  beam  waist  and  extra  space  available  in  the  long  arm.  One 
beam  waist  is  in  the  middle  of  the  crystal  (gain  medium)  and  the  other  is  at  the  center 
of  the  long  arm  and  is  available  for  elements  such  as  atomic  vapor  cells  and  nonlinear 
crystals.  For  our  specific  needs,  we  opted  for  a  ring  cavity  for  our  Ti:sapphire  laser. 

3.2  Stable  Cavity:  The  Ray  Matrix  Approach 

Once  the  choice  between  a  standing  wave  and  a  ring  cavity  has  been  made,  the 
stability  of  the  laser  cavity  must  be  addressed.  While  certain  high  power  applications 
may  require  the  use  of  an  unstable  resonator,  most  spectroscopic  applications  involve 
stable  laser  cavities.  We  address  only  this  type  of  cavities  in  this  section.  .Although 
we  look  at  the  specifics  of  a  ring  cavity,  the  general  procedure  outlined  applies  to 
standing  wave  cavities  as  well.  To  determine  the  stability  of  a  cavity,  we  must  first 
review  some  basics  from  geometrical  optics. 

Light  propagation  can  be  described  in  terms  of  rays,  which  are  geometric  curves 
representing  the  direction  of  propagation  of  energy.  In  most  cases  we  are  interested 
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CRYSTAL 


Figure  3.1:  A  basic  ring  laser  design  consisting  of  a  crystal  of  length  2f.  two  curved 
mirrors  with  radius  of  curvature  /?,  and  two  flat  mirrors  to  steer  the  beam  in  the  long 
arm.  The  distance  between  the  curved  mirrors  and  the  faces  of  the  crystal  is  d  and 
the  length  of  the  long  arm  from  Ml  to  M2  is  L.  The  crystal  has  Brewster  surfaces  at 
both  ends,  and  the  long  and  short  arms  are  taken  to  be  parallel.  .Angle  0  is  one  half 
the  fold  angle  of  the  laser  and  represents  the  angle  of  incidence  for  all  cavity  mirrors. 

in  light  propagating  predominantly  in  a  single  direction.  We  chose  this  direction  to 
be  the  ;-axis  and  refer  to  it  as  the  optic  or  cavity  axis.  .At  any  given  point  in  space, 
we  can  then  describe  a  light  ray  with  two  parameters.  The  lateral  displacement  off 
the  2:-axis  at  a  specific  ;  location  is  given  by  r(.i)  and  the  slope  of  the  ray  at  ;  is 
drfdz  =  r'(z).  If  we  form  a  column  vector  from  these  two  parameters,  we  can  obtain 
the  ray  at  a  distance  z  from  the  original  ray  (at  ^o)  by  the  transformation 


■  r(r)  ■ 

’  A 

B  ' 

r{Zo) 

.  . 

C 

D 

where  the  ABC D  transformation  matrix  (or  the  ray  matrix)  depends  on  the  optical 
element  through  which  the  ray  passes.  The  specific  parameters  for  various  optical 
elements  can  be  derived  from  the  laws  of  geometric  optics  and  are  available  in  many 
optics  and  laser  textbooks  [16.  17.  18].  The  matrices  most  often  used  in  designing 
laser  cavities  are  reproduced  in  Table  3.1  for  reference. 

The  round  trip  matrix  represents  the  transformation  of  a  ray  after  one  round  trip 
through  the  cavity.  This  matrix  is  found  by  multiplying  in  sequence  the  matrices  of 
each  optical  element  that  is  encountered  by  a  ray  in  one  round  trip.  An  analysis  of  this 
matrix  gives  us  the  stability  region,  beam  spot  size,  and  beam  waist  locations  inside 
the  cavity.  For  a  basic  ring  laser  as  depicted  in  Fig.  3.1.  we  have  two  beam  waists, 
one  at  the  center  of  the  crystal  and  one  in  the  center  of  the  long  arm.  .A  knowledge 
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Free  space  propagation  of  ciistance  d  inside 
a  medium  of  refractive  index  n. 

1  d/n 

0  1 

Thin  lens  with  focal  length  /. 

I  0 

-1//  1 

J 

Curved  mirror  with  radius  of  curvature  R 
at  an  angle  of  incidence  9,  tangential  plane. 

1  0 
—2/Rcos9  1 

Curved  mirror  with  radius  of  curvature  R 
at  an  angle  of  incidence  0.  aagittal  plane. 

1  0  ■ 
-2cosO/R  1 

Flat  interface  between  air  and  medium 
with  refractive  index  n  at  Brewster's  an¬ 
gle.  tangential  plane. 

n  0 

0  1/n 

Flat  interface  between  air  and  medium 
with  refractive  index  n  at  Brewster's  an¬ 
gle.  sagittal  plane. 

[  1  0  1 
[  0  1  J 

Table  3.1:  Transformation  Ray  Matrices  for  Optical  Elements 
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of  the  beam  waist  location  and  spot  size  is  necessary  for  mode  matching  and  efficient 
pumping  of  the  cavity  mode  volume  by  the  pump  laser.  Using  the  ray  transformation 
matrices  from  Table  3.1,  we  find  the  rovind  trip  matrix  in  the  tangential  plane  for 
starting  at  the  center  of  the  crystal  to  f»e 

\  tjn  n  0  Id  1  0 

Q  J  ~  [  0  I  J  [  0  l/«  j  [  0  1  j  [  -2IRCOS0  1 

’  1  Z.  1  f  1  0  1  r  1  d  1  r  l/n  0  1  r  1  t/n  ' 

X  '  '  .  (3.2) 

0  1  J  [  -2IRcosd  lj[0lj[0  ttj[0  1 

Here  2t  is  the  length  of  the  crystal  along  the  optic  axis,  d  is  tiie  distance  between 
either  curved  mirror  and  the  face  of  the  crystal,  and  L  is  the  total  length  of  the  long 
arm  (free  space  propagation).  For  notational  purposes,  the  subscript  t  refers  to  the 
tangential  or  the  horizontal  plane  containing  the  optic  axis  and  the  subscript  refers 
to  the  sagittal  or  vertical  plane  containing  the  optic  axis.  VVe  will  use  the  superscript 
c  to  refer  to  the  beam  waist  at  the  center  of  the  crystal  and  the  superscript  I  to  refer 
to  the  beam  waist  at  the  center  of  the  long  arm.  If  we  define  the  optical  path  length 
between  the  two  curved  mirrors  in  the  tangential  plane  as 

2lt  =  ^  +  2d.  (3.3) 

then  the  round  trip  matrix  is  simplified  and  its  elements  can  be  written  in  the  form 


For  th«'  sagittal  plane,  each  RcosO  term  is  replaced  by  RfcosO  and  the  Brewster’s 
interface  matrices  become  unit  matrices.  Defining  the  optical  path  length  Ix'tween 
the  two  curved  mirrors  in  the  sagittal  plane  as 


(3.S) 
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the  elements  of  the  round  trip  matrix  for  the  sagittal  plane  can  be  written 


■la 

fia 


2 1,  cost)  \ 


2l,cos6\ 


21,  +  L 


2l,ros0  \ 


icos0 

R 


Lco.^0\ 

«  /  ■ 


D^.  = 


id.iU) 

(d.ii) 

(d.i2) 


Once  each  element  of  the  round  trip  matrix  is  known  for  a  given  plane,  the  stability 
conditions  can  be  derived.  Since  the  ray  is  transformed  by  the  ABC D  matrix  in  each 
round  trip  through  the  cavity,  we  are  interested  in  the  range  of  the  parameters  that 
allow  rf(zo)  and  r'j{zo)  to  remain  finite  as  N  goes  to  oc  in  the  matrix  ecjuation 


v 


’  .4 

B  ‘ 

n(^o) 

. . 

C 

D 

.  . 

For  this  to  be  the  case,  the  eigenvalues  of  the  round  trip  matrix  must  be  less  than  or 
equal  to  unity.  This  leads  to  the  range  of  parameters  or  region  of  .-yt ability  defined  by 
the  inequality 

-  1  <  m  <  1  ,  (3.U) 

where  the  parameter  m  is  given  by; 


.\  complete  derivation  of  the  region  of  stability  using  eigenrays  and  eigenvalues  of  the 
round  trip  matrix  can  be  found  in  [17,  pages  599-602].  For  our  round  trip  matrix. 
.4  =  D,  so  m  =  and  the  stability  condition  becomes 


-  1  <  .4  <  1  .  (3.16) 

This  condition  must  be  satisfied  simultaneously  for  both  the  tangential  and  sagittal 
planes.  We  are  interested  in  near  confocal  spacing  between  the  curved  mirrors.  In 
order  to  investigate  the  sensitivity  of  cavity  stability  in  the  vicinity  of  confocal  spac  ing. 
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ue  iiitrodiK'*'  a  sensitivity  variable  S  by  writing  the  actual  distance  between  the  curved 
mirrors  as 

21  =  2d+  >t  =  R  +  6  .  (3.17) 


where  /?  is  the  radius  of  curvature  of  the  curved  mirrors.  In  terms  of  the  optical  pat  h 
lengths  24  and  2/,  introduced  earlier,  we  can  write 


2/f  =  R  +  6t-2t(^\- =  R  +  6t-  2ft  .  (3.18) 

and 

2l,  =  R  +  b,-  2t  (^l-^^=R  +  6,-  2fs  .  (3.1')i 

w'here  St  represents  the  6  variable  in  the  tangential  plane.  S,  represents  the  d  variable 
in  the  sagittal  plane.  We  also  introduce  the  quantities  ft  and  t,  by 


and 


(3.20) 

(3.21) 


Using  Eqs.  (3.18)  and  (3.19)  in  Eqs.  (3.4)  and  (3.9).  we  obtain  the  region  of  stability 
defined  by  Eq.  (3.16)  in  terms  of  6  for  each  of  the  planes 


2tt-  R(l  -  cose)  <  St 


<  2tt  + 

<  + 


R^cosO  —  LR[l  —  COS0) 
L  —  RcosO 
LR[  \  —  cosO) 
Lease  —  R 


(3.22) 

(3.23) 


Figure  3.2  shows  the  region  of  stability  for  both  the  tangential  and  sagittal  planes  as 
a  function  of  the  variables  S  and  L.  the  free  space  propagation  length.  The  angle  of 
incidence  is  0  =  1.5'^  and  the  radius  of  curvature  of  the  cavity  mirrors  is  /?  =  10  cm. 
Figure  3.3  shows  the  region  of  stability  as  a  function  of  the  angle  of  incidence  0 
and  the  sensitivity  variable  S  for  a  free  space  propagati  'ii  length  L  =  100  cm  and 
R  =  10  cm.  For  both  figures,  the  upper  boundary  of  stability  in  the  tangential  plane 
is  represented  by  the  dashed  line  given  by  the  equation 


St  =  2tt  + 

^maz  f  ' 


R^cose  —  LR{  \  —  COS0) 
L  —  RcosO 


(3.24) 


The  lower  bouiiflary  of  stability  in  the  tangential  plane  is  shown  by  the  clotted  line 
given  by  the  ecjuation 

=  >t,  -R{]  -co^^O)  .  i:U.j) 


The  region  between  these  two  curves  is  the  region  of  cavity  stability  for  the  tangential 
plane.  For  the  sagittal  plane,  the  upper  boundary  of  stability  is  represented  by  the 
solid  line  given  by  the  equation 


—  + 


+  L  R(  I  —  COS0) 


LcosO  —  R 

The  lower  boundary  of  the  sagittal  plane  is  given  by 

1  —  cosB  \ 


=  'it.  -  R 


cosO 


(3.26) 


(3.27) 


and  is  shown  by  the  dot-dashed  line.  The  region  between  these  two  curves  is  the 
region  of  stability  for  the  sagittal  plane. 

A  cavity  is  stable  in  both  planes  where  (he  tangential  and  sagittal  regions  of 
stability  overlap.  By  choosing  a  value  for  S  from  this  region,  a  stable  cavity  can 
be  designed  for  a  given  angle  of  incidence  9  and  free  space  propagation  length  L 
by  setting  the  distance  between  the  two  curved  mirrors  equal  to  2/  as  defined  by 
Eq.  (3.17). 

In  order  to  determine  the  spot  size  at  the  beam  waist  in  the  long  arm.  we  must 
determine  the  round  trip  matri.x  from  the  center  of  the  long  arm.  Using  a  procedure 
similar  to  that  used  to  obtain  Eqs.  (3.4)-(3.7)  and  Eqs.  (3.9)-(3.T2).  we  find  the 
elements  of  the  round  trip  matrix  starting  at  the  center  of  the  long  arm  for  the 
tangential  plane  to  be 


(3.28) 


(3.29) 


(3.30) 

(3.31) 
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Figure  3.2:  The  stability  region  for  a  ring  laser  as  a  function  of  the  sensitivity  variabl< 
6  and  the  free  space  propagation  length  L.  The  stability  region  in  the  sagittal  plane  is 
the  area  between  the  solid  and  dot-dashed  curves.  These  curves  represent  Eqs.  (3.26) 
and  (3.27)  respectively.  For  the  tangential  plane,  the  region  of  stability  is  the  area 
between  the  dashed  and  dotted  curves,  representing  Eqs.  (3.24)  and  Eq.  (3.2o). 
cavity  is  stable  in  the  region  where  the  stability  regions  of  both  planes  overlap.  For 


ANGLE  OF  INCIDENCE  (DEGREES) 

Figure  .5.3:  1  he  stability  region  for  a  ring  laser  as  a  function  of  the  sensitivity  variable 
d  and  the  angle  of  incidence  6.  The  stability  region  in  the  sagittal  plane  is  the  area 
between  the  solid  and  dot-dashed  curves.  These  curves  represent  Eqs.  (3.26)  and 
(3.27)  respectively.  For  the  tangential  plane,  the  region  of  stability  is  the  area  between 
the  dashed  and  dotted  curves,  representing  Eqs.  (3.24)  and  Eq.  (3.25).  .\  cavity  is 
stable  in  the  region  where  the  stabilities  of  both  planes  overlap.  For  these  curves,  the 
length  of  the  long  arm  Z,  is  100  cm. 
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Sitnilaily.  the  matrix  elements  in  the  sagittal  plane  are 


.4  =  1- 


= 


=  - 
D[  =  .41 


Lcos0\ 

~) 

LcosO\ 

, 

icosd 


1  - 


Al,,cosO 

ir~ 


LcosO 

R 


/.+2/J1- 


LcosO^ 

R 


R 


1 


‘2l,cos9^ 
R  . 


(3.42) 

(3.33) 

(3.34) 

(3.35) 


The  expressions  for  the  round  trip  ray  transformation  matrices  derived  in  this 
section  will  allow  us  to  determine  the  size  of  the  Gaussian  beam  waists  inside  the 
cavity.  This  information  will  be  used  to  analyze  and  compensate  for  astigmatism  in 
the  beam  profile,  as  we  now  show. 


3.3  Intracavity  Beam  Waists:  Astigmatic  Compensation 


As  seen  from  Table  3.1,  the  transformation  of  a  Gaussian  beam  at  a  curved  mirror 
and  at  an  interface  between  media  of  different  indices  of  refraction  is  different  for  the 
tangential  and  sagittal  planes.  The  tangential  plane  is  the  plane  of  incidence.  The 
sagittal  plane  is  perpendicular  to  the  plane  of  incidence.  This  difference  in  trans¬ 
formation  between  planes  is  one  source  of  astigmatism  in  laser  systems.  Nonnormal 
incidence  at  curved  mirrors  is  another  .source  of  astigmatism.  The  result  is  an  asym¬ 
metric  transverse  beam  profile,  which  is  usually  undesirable. 

In  this  section  we  describe  a  method  for  determining  the  spot  size  and  the  loca¬ 
tion  of  the  cavity  beam  waists.  We  then  describe  a  technique  for  compensating  the 
astigmatism  generated  in  the  intracavity  beam  by  the  Brewster  windows  and  curved 
mirrors  for  the  cavity  shown  in  Fig.  3.1.  .A  method  for  compensating  astigmatism  in 
a  standing  wave  cavity  composed  of  three  curved  mirrors  can  be  found  in  [19]. 

In  Section  2.4  we  saw  that  the  coefficients  of  a  ray  matrix  can  be  used  to  determine 
the  transformation  of  the  complex  Gaussian  beam  parameter  q(:).  For  a  stable  cavity, 
after  one  round  trip,  we  require  qf  =  q,  =  qo  in  Eq.  (2.61).  This  leads  to 

Aq.,  -I-  B 


= 


(  qo  A  D 


(3.36) 
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or 


Cq:  +  (D-  A)<io~  B  =  0. 


(;3.37) 


Solving  tor  l/q^  we  find 

1  D-A  ,  s/D^  -  2  AD  +  .4-=  +  \C  B 

—  =  - ±  - - - .  (3.38) 

q„  2B  2B 

Since  the  determinant  of  each  individual  ray  transfer  matrix  is  equal  to  one.  so  also 
must  be  the  determinant  of  the  round  trip  matrix 

AD-CB  =  \.  (3.:l!)) 


Using  Eq.  (3.39)  in  Eq.  (3.38)  we  find  that  the  complex  beam  parameter  q^  is  gi\en 

by  _ 


1  D-A  ,  Ji-(A  +  D)^ 
^  23  ^ '  2B 


(3.40) 


RecalJjr.g  that  1/qo  is  pure  imaginary  at  a  beam  waist,  we  see  that  A  =  D  for  the 
round  trip  matrix  starting  at  the  location  of  one  of  the  beam  waists.  For  asymmetric 
avities.  one  can  solve  for  the  location  of  the  beam  waists  by  allowing  the  starting 
point  for  the  round  trip  matrix  to  be  represented  by  a  variable  and  then  requiring 
that  ,4  =  D.  In  the  case  of  symmetric  cavities,  the  waist  location  can  be  found  by 
symmetry.  For  the  cavity  shown  in  Fig.  3.1,  the  beam  waists  are  located  at  the  center 
of  the  crystal  and  the  long  arm  as  confirmed  by  Eqs.  (3.7).  (3.12).  (3.31).  and  (3.3.3). 

.Another  method  for  defining  stable  resonators  is  suggested  by  Eq.  (3.40).  By 
requiring  that  the  beam  waist  remain  real,  we  obtain  the  condition 


{AA  D)^  <  4  . 


(3.41) 


for  stability.  This  condition  is  equivalent  to 

4  +  D 

-  1  <  —  <I  ==j>  -l<m<l,  (3.12) 

in  agreement  with  the  stability  requirement  Eq.  (3.14)  found  in  Section  3.2  using 
eigenray  analysis. 


lo  determine  the  spot  size  at  the  beam  waist,  we  recall  Eq.  (2. 12)  from  Section  2.d. 
We  then  find 

i. 


KIC 


A  Jl-{A  +  Dy 

—  —  I- - 

28 


.2 


(d.id) 


Using  the  fact  that  A  =  D.  the  spot  size  of  the  beam  at  its  waist  in  terms  of  the  ray 
matrix  elements  can  be  written  as 

in 


ii'n  = 


Wy/i  —  .4^ 


(d.ui 


Using  the  expressions  for  .4  and  8  given  in  Eqs.  (3.4).  (3.5).  (3.9)  and  (3.10).  the 
spot  size  at  the  beam  waist  at  the  center  of  the  crystal  in  the  tangential  plane  is 
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An^ 

27r 
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( —  1^  b/,  +  l(i  — 

V  RcotO  ^  ^  \l  Rco.,e) 


RcosO 


(-L-  _  i) 

V  Rcose  '■ ) 
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and  for  the  sagittal  plane  it  is 
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(  21.C036  1 

^  R  ^ 

2/,  +  L 

(1  _  31.^)' 
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27r^ 
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LcosO  1  ^ 

R  V 

cosO 

(3.15) 


(3.46) 


Similar  calculations  at  the  beam  waist  in  the  long  arm  using  the  expressions  for  .4 
and  8  given  in  Eqs.  (3.28).  (3.29).  (3.32)  and  (3.33)  yield 


(ficosS  0  Rcosei)] 

RcosO 

2tt  y 

(  Rcos9  ^  ) 

L  > 

for  the  spot  size  in  the  tangential  plane  and 
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{  Lcosff  1  ' 

1  R  h 
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^1  LcosO' 
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'2L^  -  1) 

cosO 

(3.48) 


for  the  spot  size  in  the  sagittal  plane. 

In  general,  one  finds  that  the  beam  waist  spot  sizes  in  the  two  planes  are  not 
equal  ^  ay,.  However,  we  note  that  the  spot  size  at  the  beam  waists  depend 
on  the  adjustable  variables  0.  L.  and  b.  Therefore,  by  adjusting  the  curved  mirror 


sparing  2/.  the  angle  0.  and  the  cavity  free  space  length  L.  we  can  find  a  regime  where 
astigmatism  is  minimized,  or  possibly  removed,  at  least  at  one  of  the  beam  waists.  In 
the  ideal  situation,  we  want  a  symmetric  beam  both  inside  the  crystal  and  in  the  long 
arm.  But  since  this  may  not  be  possible,  or  even  desirable  for  all  applications,  one 
must  find  the  parameters  that  best  fit  a  particular  design.  In  most  cases,  the  output 
beam  is  required  to  have  a  cylindrically  symmetric  transverse  profile.  If  the  output 
comes  from  one  of  the  flat  mirrors  steering  the  long  arm.  then  we  try  to  achieve 
it'lf  =  u'l,  in  the  long  arm  by  choosing  the  proper  values  for  6,  angle  0,  and  free  space 
length  Z,,  for  the  given  values  of  crystal  length  2t.  inde.x  of  refraction  n.  and  radius 
of  curvature  R  of  the  curved  mirrors. 

A  plot  of  the  spot  size  (beam  radius)  at  the  center  of  the  crystal  and  at  the  center 
of  the  long  arm  is  shown  in  Fig.  3.4.  The  variation  of  the  spot  size  is  shown  as  a 
function  of  the  free  space  length  L  for  an  angle  of  incidence  B  =  1.5.5'^.  Figure  3.5 
shows  the  spot  size  variation  as  a  function  of  angle  of  incidence  6  for  a  free  space 
length  L  =  100  cm.  For  both  figures,  the  radius  of  curvature  Z?  =  10  cm.  h  =  1.5  cm. 
and  A  =  780  nm,  representing  the  average  wavelength  of  la.ser  radiation  inside  the 
cavity.  The  solid  curve  shows  the  spot  size  in  the  sagittal  plane  at  the  center  of  the 
long  arm  as  given  by  Eq.  (3.48).  The  dotted  curve  shows  the  spot  size  in  the  tangential 
plane  in  the  long  arm  as  given  by  Eq.  (3.47).  For  the  beam  waist  at  the  center  of  the 
crystal,  the  spot  size  in  the  tangential  plane  is  depicted  by  the  dashed  curve  and  in 
the  sagittal  plane  by  the  dot-dashed  curve.  These  curves  represent  Eqs.  (3.45)  and 
(3.46)  respectively.  From  Figs.  3.4  and  3.5,  we  see  that  for  a  free  space  length  of 
=  100  cm.  angle  of  incidence  B  =  15..5‘^  and  b  =  1.5  cm.  the  spot  size  in  each  plane 
at  the  beam  waist  in  the  long  arm  are  equal.  Therefore,  under  these  conditions,  the 
astigmatism  in  the  long  arm  is  fully  compensated  and  the  beam  in  the  long  arm  is 
circularly  symmetric. 

Since  we  chose  to  compensate  lor  the  astigmatism  in  the  long  arm.  the  transverse 
profile  of  the  beam  inside  the  crystal  remains  asymmetric.  This  astigmatism  must  be 
kept  in  mind  when  addressing  the  problem  of  matching  the  pump  beam  mode  into 
the  cavity  mode.  This  is  the  subject  of  the  next  section. 
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Figure  3.4:  Variation  of  the  spot  size  ic^  at  the  two  intracavity  beam  waists  as  a 
function  of  L,  tfe  length  of  the  long  arm.  For  the  beam  waist  in  the  long  arm.  the 
dotted  curve  represents  in  the  tangential  plane  as  given  by  Eq.  (3.47)  and  the  solid 
curve  shows  Wo  in  the  sagittal  plane  given  by  Eq.  (3.48).  For  the  beam  waist  at  tlie 
center  of  the  crystal,  the  dashed  curve  represents  w.^,  in  the  tangential  plane  as  given 
in  Eq.  (3.4o)  and  the  dot-dashed  curve  shows  in  the  sagittal  plane  according  to 
Eq.  (3.46).  For  this  plot,  the  angle  of  incidence  is  15“.  R  =  10  cm.  6  =  1.5  cm. 
and  .\  =  780  nm. 


Figure  3.5:  Variation  of  the  spot  size  at  the  two  intracavity  beam  waists  as  a 
function  of  the  angle  of  incidence  0  for  a  free  space  propagation  length  L  =  100  cm. 
For  the  beam  waist  in  the  long  arm,  the  dotted  curve  represents  a\-,  in  the  tangential 
plane  as  given  by  Eq.  (3.47)  and  the  solid  curve  shows  u'o  in  the  sagittal  plane  given 
by  Eq.  (3.48).  For  the  beam  waist  at  the  center  of  the  crystal,  the  dashed  curve 
represents  w.^  in  the  tangential  plane  as  given  in  Eq.  (3.45)  and  the  dot-dashed  curve 
shows  Wj  in  the  sagittal  plane  according  to  Eq.  (3.46).  For  this  plot.  /?  =  10  cm. 
^  =  1.5  cm.  and  A  =  780  nm. 


29 


3.4  Mode  Matching 


Having  doterniiiipd  the  beam  waists  and  spot  sizes,  we  now  address  the  problem  (d’ 
coupling  the  pump  energy  as  efficiently  as  possible  into  the  fundamental  mode  inside 
the  cavity.  For  this  to  occur,  we  require  the  pump  beam  spot  size  and  waist  to 
coincide  with  the  cavity  mode  beam  waist  and  spot  size  as  much  as  possible.  Our 
initial  treatment  of  the  problem  assumes  a  symmetric  spot  size  at  the  beam  waist. 
In  Section  3.3  we  saw  that  this  is  not  th<-  case  for  our  cavity.  We  therefore  conclude 
this  section  with  a  technique  to  introduce  astigmatism  in  the  pump  beam  to  better 
match  the  spot  size  inside  the  crvstal.  Our  objective  is  to  optimize  the  spatial  mode 
matching  of  the  pump  and  intracavity  beams.  A  detailed  investigation  of  the  influence 
of  this  spatial  mode  matching  on  the  efficiency  of  longitudinally  pumped  solid  state 
lasers  can  be  found  in  [20]. 

The  theory  of  mode  matching  two  Gaussian  beams  is  developed  in  [loj.  The  basic- 
problem  consists  of  placing  a  focusing  lens  in  the  path  of  a  Gaussian  beam  to  generate.' 
a  second  beam  waist  at  a  predetermined  point  and  with  a  predetermined  spot  size. 
We  call  the  distance  from  the  first  beam  waist  to  the  lens  dj  and  the  distance  to 
the  second  beam  waist  di-  Using  the  transformation  matrices  from  Table  3.1.  we 
determine  the  transformation  matrix  for  a  Gaussian  beam  traveling  a  distance  of  dj. 
then  passing  through  a  lens  of  focal  length  /,  and  finally  through  a  distance  of  d2. 
The  resulting  matrix  is 


A  B  I  —  d-ij  f  di  +  d'2  —  didif  f 

C  Z)  J  ~  [  -i/f  \-d,lf 


Recalling  Eq.  (2.61).  we  relate  the  initial  beam  parameter  qi  to  the  final  beam  pa¬ 
rameter  q2  by 

_  ( 1  ~  (Izj f)  +  ^2  —  d\d2l  f  . 

Requiring  the  initial  and  final  points  of  the  transformation  to  be  beam  waists,  we  use 
Eq.  (2.43).  which  de.scribes  the  complex  beam  parameter  at  the  beam  waist,  to 
write 
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He[>‘  'i'l  rhp  spot  size  at  the  original  beam  waist  of  the  pump  beam  (inside  the 
pump  laser)  and  tcj  is  the  "transformed"  spot  size  of  the  pump  beam  at  the  second 
beam  waist  (inside  the  cavity).  Substituting  the  e.xpressions  from  Eq.  ('3.511  into 
Eq.  (3.50),  we  obtain 


■tCifCiV  1 


A 


tT{( 


,2 


A 


Equating  the  imaginary  parts  of  Eq.  (3.52),  we  find 


f-dr 


/  —  dz  u  2 
and  equating  the  real  parts  of  Eq.  (3.52)  gives 

Here  we  have  defined  the  characteristic  focal  length 

^  7rU.’iU'2 

Jo  —  ;  • 


‘-t 


(3.5:. 


3.53) 


(3.54) 


(3.55) 


To  satisfy  both  Eqs.  (3.53)  and  (3.54)  simultaneously,  we  require 


f>fo-  (3.56) 

This  means  that  if  a  pump  beam  having  a  spot  size  of  uq  at  its  original  beam  waist  is 
to  be  matched  onto  a  second  beam  waist  of  spot  size  uq,  it  must  be  focused  through 
a  lens  (or  focusing  system)  with  a  focal  length  greater  than  the  characteristic  focal 
length  defined  by  Eq.  (3.55).  VV’e  now  look  at  the  mode  matching  requirements  of  our 
laser  system. 

The  problem  for  our  laser  system  is  to  find  the  optimum  position  and  focal  length 
of  a  lens  to  focus  the  pump  beam  into  the  beam  waist  at  the  center  of  the  crystal. 
Collinear  pumping  of  the  crystal  is  desired  to  optimize  energy  transfer  from  the  pump 
source  and  simplicity  of  design.  This  requires  the  pump  to  be  focused  through  one  of 
the  cavity  curved  mirrors.  schematic  of  the  various  elements  and  distances  involved 
in  our  mode  matching  problem  are  shown  in  Fig.  3.6.  The  focusing  lens  has  a  focal 
length  /i  and  the  mirror  has  a  focal  length  of  f2-  They  are  separated  by  a  distance  a. 
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Figure  3.6:  A  schematic  of  the  optical  elements  used  in  matching  the  spot  size  of  the 
pump  beam  to  the  cavity  spot  size  at  the  beam  waist  at  the  center  of  the  crystal. 
The  lens  has  a  focal  length  fi  and  the  mirror  has  a  focal  length  of  These  elements 
are  separated  by  the  distance  a.  The  effective  lens  is  depicted  with  dash  lines  and 
has  a  focal  length  represented  by  D.  The  distance  between  the  cavity  mirror  and  the 
center  of  the  crystal  is  s  and  A  is  the  distance  between  the  effective  focusing  lens  and 
the  cavity  mirror. 


From  basic  optics  [21],  we  know  the  effective  focal  length  of  two  optical  elements  with 
focal  lengths  and  /j  separated  by  a  distance  a  is  given  by 


J _ j_  I _ ^ 

feff  fl  /2  /1/2 


(3.57) 


The  location  of  the  "effective  lens"  is  a  distance 


A  = 


q/2 

fl  +  h  ~  ^ 


(3.58) 


from  the  second  element  (cavity  mirror).  For  our  particular  application,  we  can  only 
vary  the  position  of  the  lens,  so  we  solve  Eq.  (3.57)  for 


fi  = 


f2- a 

h/ffff  -  1  ' 


(3.59) 


to  estimate  the  focal  length  requirement  of  the  matching  lens.  To  match  a  beam  waist 
with  spot  size  Wi  to  a  beam  waist  with  spot  size  W2  using  a  lens  as  depicted  in  Fig.  3.6. 
the  effective  focal  length  obtained  from  the  "focusing  system"  must  be  greater  than 
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tlu'  ( iiaracteristic  length  defined  by  Eq.  (3.55).  Because  of  space  limitations,  one 
>hould  >elect  a  lens  with  as  small  a  focal  length  as  possible  while  still  meeting  all 
other  design  requirements.  The  stability  requirements  outlined  in  Section  3.2  often 
place  stringent  requirements  on  this  focusing  lens,  since  the  cavity  mirror  must  l)e  a 
minimum  distance  from  the  center  of  the  crystal.  .Also,  very  small  changes  in  the  d 
variable  defined  in  Eq.  (3.17)  can  result  in  large  changes  in  the  spot  size  at  the  center 
of  the  crystal.  These  changes  in  turn  result  in  large  changes  in  the  mode  matching 
requirements. 

Taking  our  cavity  design  as  an  e.xarnple.  we  obtain  a  value  for  the  characteristic 
focal  length  as  follows.  Pumping  with  all  lines  from  an  argon  ion  laser,  we  take  the 
average  wavelength  to  be  A  =  500  nm.  Matching  a  beam  waist  spot  size  of  1 .8  mm 
from  the  pump  laser  to  an  average  spot  size  at  the  beam  waist  inside  the  crystal  of 
25  microns,  we  find 

%  28  cm  .  (  3.60) 

Solving  Eq.  (3.59)  assuming  =  fo  and  /2  =  -5  cm,  from  a  10  cm  concave  cavity 
end  mirror,  we  obtain  the  minimum  focal  length  for  the  lens  to  be 

«  4  +  ra  .  (3.61 ) 

VV  e  see  that  the  minimum  focal  length  depends  on  the  minimum  acceptable  spacing 
between  the  lens  and  the  cavity  mirror.  .Assuming  we  want  a  minimum  of  5  cm  from 
the  mirror  to  the  lens  for  mounting  and  adjusting,  we  have  a  minimum  focal  length 
of 

f,  ss  8  cm  (3.62) 

for  the  focusing  lens  in  our  mode  matching  problem.  For  our  laser  system,  a  10  cm 
achromatic  lens  was  selected  to  meet  the  minimum  fi  requirements,  and  allow  for 
variations  in  the  pump  laser,  the  value  of  S.  spot  sizes  and  cavity  configurations. 

The  position  of  the  "effective  lens"  from  the  beam  waist  inside  the  crystal  is 
determined  from  the  formula 

D  =  A  +  »  =  (M:!) 
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where  ir^  is  being  matched  to  it’2-  The  distance  .s  between  the  cavity  mirror  and  the 
center  of  the  crystal  is  determined  by  the  value  of  6  according  to  Eq.  (3.17).  Sirue 
this  value  is  determined  from  the  stability  requirements  of  Eqs.  (3.2'2)  and  (3.23).  we 
consider  it  to  be  a  constant  for  the  mode  matching  problem.  The  value  of  s  is  given 
by 

s  =  ^(R  +  6).  (3.61) 

where  R  is  the  radius  of  curvature  of  the  cavity  mirror  in  the  focusing  system  and 
is  chosen  from  the  stability  diagrams  (see  Figs.  3.2  and  3  3). 

To  find  the  position  of  the  focusing  lens,  we  first  solve  Eq.  (3..")7)  for  /e//  explicitly 


Jeff  = 


^1/2 

fi  +  h  —  (I 


(3.6.)) 


Using  Eqs.  (3.65)  and  (3.58)  we  rewrite  Eq.  (3.63)  as 


a/2 


/i  +  /2  —  a  f\  f 2  —  a  fi  —  a 


fifz 


W2 


/1/2 


-  P 

J  0 


(3.66) 


In  most  applications,  u’2  -C  tt'i  and  /e//  «  /o.  so  we  neglect  the  second  term  and 
solve  for  a 

^  flh  —  +  fz)  /1/2  —  2(^  T  +  .^2) 

- —s - =  - iPW^) - ■ 

For  a  10  cm  lens  (/i),  a  curved  mirror  of  R  =  10  cm  (fi  =  —5  cm)  and  S  =  1.5  cm 
we  obtain 

a  %  7.3  cm  .  (3.68) 


Thus  by  placing  the  10  cm  lens  a  distance  of  7.3  cm  behind  the  cavity  end  mirror,  we 
can  focus  the  pump  beam  to  its  minimum  spot  size  at  the  center  of  the  crystal.  In 
pra-'tice.  this  calculated  distance  is  a  “starting  point"  for  the  laser  design  and  shoidd 
be  adjusted  by  trial  and  error  to  compensate  for  the  term  neglected  in  Eq.  (3.66)  and 
optimize  the  performance  of  the  laser. 

We  saw  from  Section  3.3  that  the  spot  size  at  the  center  of  the  crystal  may  be 
asymmetric.  To  effectively  mode  match  the  pump  volume  and  the  intracavity  volume 
inside  the  crystal,  we  need  to  introduce  astigmatism  into  the  pump  beam.  This  can  be 
done  by  rotating  the  focusing  lens  about  the  vertical  axis  by  a  few  degrees.  In  fact. 


some  astigmatism  is  already  present  due  to  the  nonnormal  approach  of  the  pump 
!)eam  oti  the  curved  cavity  mirror.  Calculations  by  Schulz  [22]  for  a  14''  angle  of 
incidence  at  the  curved  cavity  mirror  show  a  required  rotation  of  approximately  S\ 
Our  experience  showed  that  a  similar  rotation  in  the  focusing  lens  r=‘sulted  in  increased 
output  power. 

In  this  section  we  have  described  the  problem  of  matching  the  mode  of  the  pump 
beam  to  that  of  the  cavity.  VVe  saw  that  this  can  be  accomplished  by  proper  sele<'tion 
and  placement  of  a  focusing  lens  outside  the  cavity.  \Ve  also  briefly  discussed  how 
astigmatism  may  be  introduced  in  the  pump  beam  to  better  match  the  cavity  mode 
volume.  We  now  investigate  various  intracavity  elements  that  can  be  used  in  the  laser 
system. 


3.5  Intracavity  Elements 

The  specific  application  for  a  laser  system  dictates  the  type  of  intracavity  elements 
recjuired.  In  this  section,  we  describe  some  intracavity  elements  that  can  be  used 
inside  a  laser  cavity  to  meet  design  requirements.  These  include  unidirectional  lasing 
and  single  wavelength  operation.  Our  objective  in  this  section  is  only  to  give  a  brief 
overview  of  certain  intracavity  elements.  A  more  rigorous  treatment  of  the  theory  of 
these  elements  can  be  found  for  example  in  [23]. 

Since  some  of  these  elements  rely  on  the  birefringent  characteristics  of  certain 
materials,  we  begin  this  section  with  a  basic  review  of  birefringence. 

3.5.1  Birefringence 

Birefringence,  or  double  refraction,  is  a  property  of  anisotropic  crystals  where  the 
index  of  refraction  seen  by  incident  light  is  a  function  of  its  polarization.  Refractive 
indices  are  determined  from  the  dielectric  permeability  of  the  medium.  Here  we  only 
consider  uniaxial  crystals.  In  these  crystals,  two  axes  defined  by  the  crystalline  struc¬ 
ture  have  the  same  dielectric  permeability  while  the  third  axis  has  a  different  value. 
The  axis  with  the  unique  dielectric  permeability  is  referred  to  as  the  optic  axc-i  and 
represents  a  unique  direction  of  propagation  inside  the  crystal.  The  electromagnetic 


waves  representing  the  light  incident  on  a  birefringent  medium  are  broken  down  into 
components  based  on  the  orientation  of  the  E  field  vector  relative  to  the  o|)tic  axis. 
The  components  of  the  electromagnetic  waves  with  the  E  field  vector  oriented  along 
the  optic  axis  experience  one  index  of  refraction  n,  and  travel  at  the  speed  r/n,.  The 
trajectory  of  these  waves  represents  the  fxtraordinary  btam.  The  components  of  the 
electromagnetic  waves  with  the  E  field  vector  orthogonal  to  the  optic  axis  experience 
a  different  index  of  refraction  and  travel  at  the  speed  c/n,,.  The  trajectory  of  these 
waves  represent  the  ordinary  beam. 

The  different  speeds  of  the  ordinary  and  extraordinary  beams  result  in  a  phase 
difference  between  the  two  types  of  electromagnetic  waves.  This  phase  difference  is 
represented  by 


(:Tf)9) 


where  d  represents  the  thickness  of  the  birefringent  medium.  By  varying  this  thickness 
one  can  control  the  phase  difference  induced  by  the  birefringence.  If  the  overall  phase 
difference  is  d  =  ±r/2.  then  the  birefringent  crystal  changes  linearly  polarized  light 
into  circularly  polarized  light,  provided  the  crystal  is  oriented  so  that  the  optic  axis 
is  at  45°  to  the  plane  of  polarization  of  the  incident  beam.  This  optical  element 
is  known  as  a  quarter  wave  plate.  If  the  overall  phase  difference  is  d  =  ±7r.  then 
the  birefringent  crystal  rotates  the  plane  of  polarization  by  an  angle  of  2o.  where 
O  represents  the  angle  between  the  plane  of  polarization  and  the  optic  axis  of  the 
crystal.  Such  optical  elements  are  known  as  polarization  rotators  or  half  wave  plate.'^. 
The  direction  of  rotation  is  determined  by  the  direction  of  propagation  through  the 
birefringent  crystal. 


3.5.2  Unidirectional  Lasing 

.A  basic  ring  cavity,  as  shown  in  Fig.  .3.1,  allows  for  an  electromagnetic  wave  to 
propagate  in  both  directions.  This  means  that  as  the  gain  medium  is  pump<’<l  and 
energy  inside  the  gain  medium  builds  up.  the  modes  of  each  direction  compete  against 
each  other.  This  competition  can  result  in  bistability  or  chaotic  behavior.  Many 
applications  for  ring  lasers  rerjuire  the  laser  to  be  stable  and  therefore  the  dir<‘ction 
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ot  lasing  to  be  controlled.  Thi  bi  usually  accomplished  by  giving  the  cavity  a  favored 
direction  bv  suppressing  or  enhancing  one  of  the  directional  modes. 

One  technique  for  introducing  a  favored  direction  inside  the  cavity  is  to  reflect  the 
output  from  one  direction  of  the  cavity  back  onto  itself  with  an  external  mirror.  I  bis 
method  has  the  effect  of  transferring  a  percentage  of  the  energy  from  the  undesired 
directional  mode  to  the  desired  directional  mode,  '^•te  drawback  of  this  technique 
is  that  pure  unidirectional  lasing  cannot  be  obtained  since  energy  buildup  in  the 
unfavored  direction  inside  the  cavity  is  not  fully  suppressed.  Furthermore,  this  type 
of  feedback  may  lead  to  phase  and  frequency  instability. 

To  obtain  unidirectional  lasing,  one  must  suppress  the  undesired  directional  cavitv 
mode  sufficiently  to  ensure  that  the  threshold  condition  for  that  mode  is  not  met . 
This  suppression  is  often  accomplished  with  an  optical  diode  consisting  of  a  Faraday 
rotator  and  some  type  of  optical  compensator.  The  manner  by  which  these  elements 
produce  unidirectional  lasing  is  as  follows. 

Laser  systems  usually  contain  surfaces  oriented  at  Brewster's  angle  with  respect  to 
the  cavity  axis  in  order  to  minimize  the  insertion  loss  of  intracavity  elements.  These 
Brewster  windows  give  rise  to  a  dominant  polarization  of  the  electromagnetic  wave 
traveling  inside  the  cavity.  This  dominant  polarization  is  in  the  tangential  plaiie.  or 
the  plane  of  incidence  defined  by  the  cavity  axis  and  the  normal  to  the  surface  of  the 
Brewster  window.  This  polarization  component  experiences  no  loss  at  the  Brewster 
interface. 

Certain  materials  rotate  the  plane  of  polarization  of  the  electromagnetic  wave  as 
it  passes  through  the  medium.  Birefringent  crystals,  alluded  to  in  Section  3.5.1.  are 
such  materials.  Other  materials  rotate  the  polarization  of  an  electromagnetic  wave 
when  they  are  placed  in  a  magnetic  field.  This  effect  is  known  as  the  Faraday  ffffct 
and  elements  using  this  principle  are  referred  to  as  Faraday  rotators.  The  extent  of 
rotation  of  polarization  per  unit  length,  or  the  “rotary  power,  "  of  a  Faraday  rotator 
is  proportional  to  the  intensity  of  the  magnetic  field  and  an  intrinsic  property  of 
the  material  represented  by  the  Verdet  constant.  I'nfortunately.  materials  with  large 
\erdet  constants  usually  also  have  high  absorption  coefficients.  This  means  that  large 
rotation  angles  can  only  be  obtained  at  the  expense  of  high  transmission  losses.  The 


(lirtH  tion  of  polariiiatioii  rotation  is  determined  solely  by  the  direction  of  the  magnetic 
fiehl  anrl  is  therehne  independent  of  the  direction  of  propagation.  In  other  words, 
if  the  polarizatit)!!  of  a  wave  changes  by  an  angle  o  after  passing  once  through  a 
Faraday  rotator,  it  will  have  a  total  rotation  of  'lo  if  it  is  reflected  back  through  the 
same  Faraday  rotator  in  the  opposite  direction. 

Fo  induce  loss  in  one  direction  inside  the  cavity  and  not  the  other,  a  Faraday 
rotator  is  combined  with  a  half  wave  plate.  The  half  wave  plate  is  oriented  to  com¬ 
pensate  for  the  rotation  induced  by  the  rotator.  In  one  direction,  the  half  wave  plate 
rotates  the  polarization  back  to  its  original  oiientation.  For  this  direction,  the  wave 
suffers  no  loss  as  it  travels  through  the  Brewster  windows  inside  the  cavitv.  In  the 
other  direction  however,  the  waveplate  adds  to  the  rotation  of  polarization  and  the 
electromagnetic  wave  is  attenuated  by  traveling  through  a  Brewster  window.  Esen 
with  small  rotation  angles,  the  attenuation  of  the  wave  is  usually  sufficient  to  fully 
suppress  the  undesired  directional  cavity  mode. 

3.5.3  Tuning  to  Single  Wavelength 

Lasers  operating  with  broad  gain  medium  require  some  type  of  wavelength  selecting 
device  if  single  wavelength  operation  is  desired.  To  fully  utilize  broad  t unability, 
we  desire  single  wavelength  operation  with  minimal  loss  to  the  desired  wavelength. 
VVe  also  want  to  have  a  high  transmission  bandwidth  at  the  selected  wavelength  and 
maximum  rejection  at  the  unwanted  wavelengths.  The  optical  element  best  suited 
for  these  requirements  is  the  birefringent  filter. 

Since  their  introduction  as  a  tuning  element  for  broadband  lasers,  the  theory  and 
design  of  birefringent  filters  has  been  an  area  of  active  research  [24.  25.  26.  27].  These 
filters  have  been  used  as  tuning  elements  in  a  number  of  broadband  lasers,  including 
dye.  high  power  Nd:gla.ss.  Ti:sapphire.  and  He^Se"^  [28.  29].  The  basic  theory  behind 
birefringent  filters  is  as  follows. 

.\s  was  shown  in  Section  3.5.1.  birefringent  crystals  can  be  used  to  rotate  the 
polarization  of  an  incident  electromagnetic  wave.  .A  birefringent  plate  which  can  be 
used  as  a  wavelength  selecting  devise  is  shown  in  Fig.  3.7.  The  birefringent  plate  is 
placed  into  the  laser  cavity  such  that  the  normal  of  the  plate  is  at  Brewster's  angle 
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Figure  3.7:  A  birefringent  plate  used  as  an  intracavity  tuning  element.  The  normal 
to  the  plate  is  at  Brewster's  angle  $b  to  the  cavity  axis.  The  optic  axis  of  the 
crystal  makes  the  angle  ^  with  the  path  of  the  refracted  beam  inside  the  crystal.  The 
thickness  of  the  plate  is  d. 


6b  to  the  cavity  axis.  The  refracted  beam  inside  the  plate  travels  in  a  direction  that 
makes  the  angle  ^  with  the  optic  axis  of  the  crystal.  The  phase  difference  induced 
between  the  ordinary  beam  and  the  extraordinary  beam  of  wavelength  A  after  it  has 
travelled  through  the  plate  of  thickness  d  is  given  by 


^  _  2ir{n.  —  n^)d$in^^ 

XsinOs 


(3.70) 


Electromagnetic  waves  for  which  this  phase  difference  is  a  multiple  of  27r  will  experi¬ 
ence  no  rotation  in  polarization  and  hence  no  attenuation  as  they  pass  through  the 
end  surface  of  the  plate  at  Brewster's  angle.  The  wavelength  of  these  waves  can  be 
determined  by  setting  d  =  2mir  in  Eq.  (3.70)  and  solving  for  A 


A  = 


(He  -  no)dsin^^ 


(3.7i; 


msinOB 

where  rn  is  an  integer.  Since  the  plate  is  introduced  into  the  cavity  at  Brewster's 
angle,  all  other  wavelengths  will  suffer  loss  as  they  exit  the  plate. 

In  practice,  birefringent  filters  are  composed  of  a  series  of  birefringent  plates. 
This  is  to  provide  stronger  suppression  of  the  undesired  wavelengths  and  narrowing 
the  transmission  bandwidth  of  the  filter.  By  rotating  the  filter  around  the  normal 
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to  tho  surface  of  the  plate  and  hence  keeping  the  plate  at  Brewster  s  angle  to  the 
car  ity  axis,  one  can  change  the  angle  ^  between  the  optic  axis  of  the  crystal  and  the 
direction  of  the  refracted  beam  inside  the  plate.  From  Eq.  (d.Tl ).  we  see  that  such  a 
change  in  ^  will  result  in  a  change  in  the  "privileged"  wavelength  which  is  transmittefl 
without  loss.  In  this  manner,  one  can  tun»'  the  laser  to  a  desired  wavelength. 

.A  more  rigorous  treatment  of  the  theory  and  design  of  these  filters  can  be  found 
in  the  references  cited  at  the  beginning  of  this  subsection.  For  a  study  of  cavity 
modes  containing  birefringent  filters  see  [30]  and  a  detailed  discussion  of  the  stability 
of  birefringent  plates  can  be  found  in  [31]. 

In  this  chapter,  we  have  disctissed  the  basic  principles  of  designing  a  laser  system. 
VVe  discussed  the  advantages  and  disadvantages  of  standing  wave  and  ring  cavities. 
We  presented  an  approach  to  the  cavity  design  in  terms  of  ABC D  ray  transformation 
matrices  and  showed  how  they  can  be  used  to  analyze  the  stability  of  laser  cavities. 
We  also  developed  a  technique  for  compensating  for  astigmatism  in  the  intracavity 
beam.  We  then  discussed  the  mode  matching  problem  and  showed  how  to  optimize  the 
pump  mode  coupling  into  the  cavity  mode  inside  the  gain  medium.  We  also  presented 
a  brief  survey  of  the  intracavity  elements  that  can  be  used  for  unidirectional  lasing 
and  tuning  the  laser  to  a  single  wavelength.  With  an  understanding  of  the  basics  of 
laser  design,  we  now  turn  our  attention  to  the  theory  of  laser  dynamics. 
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Chapter  4 


Laser  Dynamics 


An  understanding  of  the  behavior  of  any  dynamical  system  requires  a  theoretical 
framework  within  which  the  system  dynamics  can  be  analyzed.  Such  a  framewink 
usually  consists  of  a  small  number  of  ecjuations  of  motion.  These  equations  of  motion 
are  derived  for  an  idealized  system  that  mimics  the  physical  system  under  investiga¬ 
tion  as  much  as  possible.  In  this  chapter,  we  present  two  basic  equations  of  motion  for 
a  laser  system  known  as  the  laser  rale  equations.  VVe  then  show  how  these  equations 
can  be  used  to  analyze  relaxation  oscillations  exhibited  by  certain  lasers.  Many  solid 
state  lasers  exhibit  these  oscillations.  In  particular,  the  Ti:sapphire  laser  riescribed 
in  this  thesis  exhibits  pronounced  relaxation  oscillations  when  perturbed  by  external 
excitations  from  the  steady  state.  We  also  describe  a  numerical  method  for  modeling 
the  dynamical  behavior  of  the  laser  within  the  rate  equation  approximation.  Finally, 
we  look  at  the  threshold  condition  for  laser  action  and  see  how  this  condition  can  be 
used  to  determine  the  intrinsic  cavity  losses  of  a  laser  system. 

4.1  Laser  Rate  Equations 

The  treatment  of  laser  dynamics  using  rate  equation  theory  involves  the  exchange 
of  energy  between  two  important  dynamical  variables,  the  number  of  photons  inside 
the  laser  cavity  q  and  the  number  of  excited  atoms  n  inside  the  laser  gain  medium. 
These  two  quantities  evolve  in  time  according  to  the  coupled  cavity  and  atomic  rate 
equation-^.  The  derivation  of  these  equations  can  be  found  in  most  textbooks  on  la.ser 
theory  [16.  17.  .32.  33].  The  form  of  these  equations  depends  somewhat  on  the  model 
used  to  describe  the  gain  medium.  VVe  opt  for  the  four  level  laser  model  for  the 
Tiisapphire  gain  medium,  as  shown  in  Fig.  (4.1).  We  assume  that  the  upper  pump 
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Figure  4.1;  A  general  model  of  the  four  level  laser  pumping  scheme.  The  decay  rates 
from  the  upper  pump  band  ^3  and  the  lower  lasing  level  7,  are  assumed  to  be  much 
greater  than  the  decay  rate  of  the  upper  lasing  level  72.  Rp  represents  the  rate  at 
which  atoms  are  introduced  into  the  upper  lasing  level. 


band  and  the  lower  lasing  level  decays  instantaneously  compared  to  the  fluorescence 
lifetime.  Assuming  single  mode  laser  operation,  we  arrive  at  the  cavity  rate  equation, 
which  describes  the  time  evolution  of  the  number  of  photons  q  inside  the  laser  cavity 

=  A'n(<)[g(<)  +  1]  -  7,9(0  .  (4.1) 


Here  K  is  the  spontaneous  emission  rate  per  mode  and  7,  is  the  cavity  decay  rate 
w-hich  is  related  to  the  round  trip  cavity  loss  0  (in  percent  per  pass)  and  the  cavity 
round  trip  time  T  by 


Ic  — 


loor 


=  T. 


(4.2) 


Eq.  (4.2)  introduces  r,  as  the  photon  lifetime  inside  the  cavity.  The  “1"  inside  the 
square  brackets  in  Eq.  (4.1)  is  sometimes  referred  to  as  the  "extra  photon"  and 
represents  the  first  photon  inside  the  cavity  mode  due  to  the  spontaneous  emission 
process.  Without  this  spontaneously  emitted  photon,  the  cavity  photons  would  not 
"build  up."  Therefore,  this  “1"  is  required  when  Eq.  (4.1)  is  used  to  describe  the 
initial  build  up  of  laser  action.  .At  steady  state,  or  other  conditions  where  q  >  1.  this 
term  can  be  ignored. 
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1  h('  atomic  rate  equation,  which  describes  the  number  of  atoms  n  in  the  u,/pcr 
state  is  given  Iw 

(ln{t) 


(4.3) 


df  ■ 

where  Rp  represents  the  rate  at  which  atoms  are  introduce  into  the  upper  lasing  level 
and  -)2  is  the  upper  lasing  level  decay  rate.  The  lifetime  of  the  upper  state,  r,.  is 
related  to  its  decay  rate  by 

(i.n 


.  —  .,--1 
!2  -  '2 


Since  it  is  assumed  that  there  is  no  accumulation  of  atoms  in  the  lower  lasing  level. 
Eq.  (4.3)  also  describes  the  evolution  of  laser  inversion  or  the  population  inversion. 
The  steady  state  solution  of  Eq.  (4.1)  and  Eq.  (4.3)  gives  us  an  expression  for  the 
number  of  atoms  in  the  upper  state 


- 


fc 


l<  ' 


and  the  number  of  photons  inside  the  cavity 


—  i  1 


Ic 


Here  we  have  introduced  the  laser  pump  parameter  r  by 


r  = 


RpK 

~t2lc 


(4.0) 


(  l.h) 


(4.i 


The  laser  pump  parameter  can  also  be  w’ritten  as  the  ratio  of  the  pump  power  to  the 
threshold  pump  power 


Wp 

IV', 


th 


Note  that  at  threshold,  r  =  I  and  from  Eq.  (4.6)  we  can  write 

K  -pi 


(4.S) 


(4.9) 


Comparing  Eq.  (4.9)  with  (4.o).  we  see  that  in  the  steady  state,  la.ser  inversion  is 
"clamped"  at  its  threshold  value 


I.IO) 


We  use  these  rate  equations  to  look  at  one  specific  dynamical  behavior  known  as 
relaxation  oscillations  exhibited  bv  certain  la.sers. 


4.2  Relaxation  Oscillations:  Laser  Rate  Equation  Analysis 


M  any  dynaniir  sxstems  operating  in  the  steady  state  will  return  to  their  steady  state 
after  being  perturbed  if  the  steady  state  is  stable.  The  approach  back  to  the  steady 
state  may  not  be  monotonic.  Instead,  the  system  may  approach  the  steady  state 
in  an  oscillatory  manner.  These  oscillations  are  usually  damped  and  are  known  as 
relaxation  oscillations.  In  this  section,  we  analyze  the  rate  equations  of  .Section  4.1  in 
the  presence  of  perturbations  from  steady  state  and  see  under  what  conditions  a  las('r 
system  e.xhibits  rela.xation  oscillations.  We  will  determine  the  ecpiations  of  motion 
for  the  relaxation  oscillations  in  terms  of  characteristic  laser  parameters. 

Since  the  laser  rate  equations  (4.1)  and  (4.3)  are  coupled,  a  perturbation  in  either 
the  number  of  photons  inside  the  cavity  or  the  number  of  atoms  in  the  upper  state  will 
result  in  a  perturbation  of  the  other.  Treating  the  perturbations  as  small  compared 
to  the  steady  state  values,  we  look  for  solutions  to  Eqs.  (4.1)  and  (4.3)  of  the  form; 


<'/(0  =  <-(0  where  c(/)  <  .  (4.11) 

n(0  =  n„  +  ri{t)  where  T)(t)  <  n„  .  (4.12) 


Substituting  Eqs.  (4.11)  and  (4.12)  into  Eqs.  (4.1)  and  (4.3)  and  keeping  only  terms 
linear  in  e  and  r/.  we  find  the  coupled  rate  equations  for  the  perturbations  to  be 

de{t) 


dt 

dv(t) 

dt 


=  (r- 1)7-27(0 
=  -  nx^iit) 


(4.13) 

(4.14) 


If  we  take  the  time  derivative  of  Eq.  (4.13)  and  use  Eq.  (4.14)  for  i).  we  obtain  a 
closerl  equation  for  e(/). 


dt^ 


Mt)  , 


(4.15) 


A  'similar  ecpiation  is  obtained  for  population  perturbation.  Equation  (4.15)  describes 
the  perturbation  in  the  cavity  photon  number.  By  introducing  parameters  d  and  u,',, 
as 


23  =  r~;2  and  =  (r  —  l)"2*r  •  (4.161 
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we  can  rewrite  the  equation  of  motion  for  the  perturbation  of  cavity  photon  number 
as 

d‘f{l)  ) 

— +  2.J— +  -■„<(()  =0.  H.I7I 

This  is  the  eciuation  of  motion  for  a  damped  harmonic  oscillator  with  natural  fre¬ 
quency  and  energy  decay  constant  2T.  The  solution  of  this  equation  is  of  the 
form 

c(/)  =  6,e— ^  (i.lS) 

Substituting  Eq.  (4.18)  into  Eq.  (4.17)  we  get  the  auxiliary  equation 

•jj"  +  =  0  . 

which  determines  Solving  this  ecpiation  for  u)  we  find 

u!  =  —i,3  ±  —  3^  =  —i3  i  wvV  , 

where  we  have  introduced  a  shifted  frequency  u;^  by 


1.1b) 


(4.20) 


(4.21 : 


.As  is  well  known,  dissipation  shifts  the  frequency  of  oscillation  [34]. 

The  behavior  of  f(<)  depends  on  the  relative  magnitude  of  u,'o  and  3.  We  consider 
three  cases  separately. 

a.)  For  the  case  u;^  <  3^.  u;  is  pure  imaginary  and  the  solution  from  Eq.  (4.18) 
takes  the  form 

e(f)  =  .  (4.22) 

Thus  any  perturbation  from  the  steady  state  exponentially  decays  to  zero.  The  system 
is  said  to  be  overdamped.  In  this  case  no  oscillations  in  the  laser  photon  number  will 
be  observed.  The  condition  <  3^  can  also  be  expressed  in  terms  of  decay  rates 
of  the  la.ser  system.  By  recalling  the  definitions  of  3  and  from  Eq.  (4.16).  we  find 
that  the  condition  for  an  overdamped  response  is 


(r  -  1)0.-  < 


r^2 


(4.23) 
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riiis  implies  that  the  cavity  decay  rate  must  be  less  than  or  of  the  same  order  of 
magnitude  as  the  upper  state  atomic  decay  rate.  This  is  characteristic  of  most  gas 
lasers,  where  relaxation  oscillations  are  not  observed. 

b.)  When  =  P.  the  two  roots  to  Eq.  (4.19)  are  equal.  For  this  case,  the 
solution  from  Eq.  (4.18)  takes  the  form 

e{t)  =  (At  +  B)  .  (4.21) 


and  the  system  is  said  to  be  critically  damped.  Any  perturbation  in  the  system  decays 
to  zero  asymptotically  with  time.  This  condition  exits  in  a  laser  system  when 


(r-  l)7c  =  (4.2.-)) 

or  the  cavity  decay  rate  and  upper  state  decay  rate  are  of  the  same  order  of  magnitude. 

c.)  For  the  case^^,’^  >  P ,  is  real.  Using  Eq.  (4.20).  the  solution  from  Eq.  (4.18) 

takes  the  form 

e(t)  =  e"'^'  .  (4.26) 

This  solution  can  also  be  rewritten  in  the  form 


e(t)  =  e.^e  ^‘cos{u;rt  +  <p)  . 


(4.27) 


In  this  case,  the  system  is  said  to  be  undtrdamped  and  displays  an  exponentially 
damped  sinusoidal  response  to  the  perturbation  as  it  returns  to  steady  state.  In 
terms  of  decay  rates,  we  expect  to  see  these  relaxation  oscillations  in  a  laser  system 
where 

4c  >  72.  (4.28) 

This  is  the  case  for  many  solid  state  lasers.  Expressed  in  terms  of  laser  parameters 
from  Eq.  (4.16)  the  oscillations  around  the  steady  state  have  a  frequency  given  by 


The  perturbation  has  an  exponentially  decaying  envelop  characterized  by  a  decay  rate 
given  by 

(4.40) 


46 


From  Fqs.  (4.29)  and  (4.'i0)  vve  see  that  one  can  determine  certain  laser  parameter'^ 
by  analyzing  the  rela.xation  oscillations.  Indeed,  this  was  one  of  the  objectivf's  of  t  he 
e.xperiments  discussed  in  Chapter  5. 

Once  we  understand  the  characteristic  equations  of  motion  for  a  dynamic  system 
and  their  controlling  parameters,  we  can  now  address  the  question  of  whether  these 
equations  adequately  describe  the  system  dynamics  in  general.  This  is  the  subject  of 
the  next  section. 


4.3  Relaxation  Oscillations:  Numerical  Modeling 


One  of  the  best  ways  to  study  a  dynamic  system  is  to  model  it  and  compare  the 
observed  response  with  that  predicted  by  the  model.  In  this  section,  we  describe  a 
method  for  modeling  the  relaxation  oscillations  exhibited  by  a  laser. 

Although  the  coupled  laser  rate  equations  of  Section  4.1  are  fairly  simple  in  form, 
exact  analytical  solutions  to  them  are  limited  and  complex  in  form  [35].  Therefore, 
in  order  to  effectively  model  any  type  of  laser  dynamic  these  rate  equations  must 
be  solved  numerically.  The  number  of  photons  q  and  the  number  of  excited  atoms 
n  are  very  large  for  realistic  systems.  These  large  numbers  are  difficult  to  manage 
numerically.  We  therefore  introduce  scaled  variables  that  are  better  suited  for  nu¬ 
merical  modeling.  The  scaling  should  be  chosen  in  such  a  way  as  to  minimize  the 
number  of  parameters  in  the  equations  of  motion.  Also,  the  scaling  should  allow  the 
equations  to  be  solved  with  parameters  measured  from  the  experiment.  We  present 
two  different  procedures  for  scaling.  These  are  suitable  for  modeling  different  types 
of  perturbations. 

For  the  system  when  gain  is  perturbed  we  scale  the  cavity  photon  number  and 
the  number  of  excited  atoms  by  introducing  the  variables 


j.  ^  ±  =  111 

fl3  72 
n  K  n 

nth  7.- 


(4.31) 

(4-32) 


Here  q,  represents  the  i^atiiration  photon  number  for  the  active  atoms.  It  is  defined 
as  the  ratio  of  population  decay  rate  ^nd  the  spontaneous  emission  decay  rate  per 
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mode  l\  bv 


(l.d:5) 


The  scale  factor  tor  population  nth  is  the  thixshold  population  inversion  defined  by 
Fiq.  (4.9).  Equations  (4.31)  and  (4.32)  imply  that  x  is  the  cavity  photon  number 
in  units  of  saturation  photon  number  and  y  is  the  population  inversion  in  units  of 
threshold  inversion.  Using  Eq.  (4.7)  we  can  e.xpress  the  population  pumping  rate  in 
terms  of  laser  pump  parameter  r  and  various  decay  rates  as 


(4.3  n 


Substituting  the  scaled  variables  from  Eqs.  (4.31).  (4.32)  and  (4.34)  into  Ef{s.  (4.1) 
and  (4.3).  the  rate  equations  for  the  scaled  variables  are  found  to  be 


dx 

It 

dt 


ic  \y 


=  liir  -y  -  yx)  . 


Einally.  scaling  time  in  terms  of  the  cavity  lifetime  by 

t 

r  =  7c<  =  —  . 

r,- 


we  can  rewrite  Eqs.  (4.35)  and  (4.36)  as 


dx 

dr 

dr 


(y  -  l)x  +  — 


s(r  -  y  -  xy)  . 


(4.35) 

(4.36) 


(4.37) 


(4.38) 

(4.39) 


Here  s  is  the  ratio  of  the  population  decay  rate  72  and  the  cavity  decay  rate  7,- 


—  .  (4.40) 

Ic 


These  equations  are  appropriate  for  modeling  a  system  where  the  gain  of  the  laser 
is  modulated.  .Modulation  of  gain  can  be  described  in  terms  of  the  pump  parameter 
alone. 
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l  or  systems  where  loss  is  modulated,  both  and  r  are  changing,  so  a  different 
scaling  is  more  appropriate.  For  dynamical  stiidies  with  this  type  of  system,  we 
replace  y  dehned  iti  Eq.  (4.32)  with  the  variable 


An 


(4.11  ) 


72 

In  this  model  we  scale  both  the  number  of  atoms  in  the  upper  state  and  the  number 
of  photons  in  the  cavity  by  the  saturation  photon  number  defined  by  Eq.  1 4.33).  We 
scale  the  time  in  terms  of  the  upper  state  lifetime  by 


t 

r  =  7-2/  =  —  . 

r-i 


(4.42) 


Using  Eqs.  (4.31).  (4.41)  and  (4.42)  into  Eqs.  (4.1)  and  (4.3)  we  obtain  the  new  rate 
equations 


dx  /  1  \ 

dr 

dz 

—  =  R-z{x  +  l). 


where  we  have  introduced  a  new  parameter  p  by 


(4.43) 
I  1.11) 


The  normalized  pumping  rate 


R  = 


RpK 

7| 


=  rp 


(l.lo) 


(1.16) 


then  depends  only  on  the  pumping  rate  Rp  and  stays  constant  as  long  as  Rp  is 
(xmstant.  The  only  parameter  that  changes  when  the  cavity  loss  is  modulated  is  p. 
defined  in  Eq.  (4.45).  It  is  helpful  to  solve  Eqs.  (4.43)  and  (4.44)  for  the  steady  state. 
For  the  scaled  cavity  photon  number  at  steady  state  we  find 


S3  — 


(  1.17) 


For  the  scaled  population  inversion  we  find  (neglecting  the  \/q,  term) 


=  P  ■ 


(4.48) 
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WV  will  fiiul  these  steady  state  values  useful  in  Section  5.2.1  when  we  use  Eqs.  i  l.l.lj 
and  (  1.14)  to  model  our  experimental  observations. 

I  sing  Eqs.  (4.58)  and  (4.59).  or  Eqs.  (4.45)  and  (4.44).  one  could  model  a  system 
numerically  depending  on  whether  the  gain  or  the  loss  is  modulated.  These  e(niations 
are  scaled  in  such  a  manner  that  only  a  single  parameter  is  needed  to  describe  the 
modulation  of  the  system.  For  Eqs.  (4.38)  and  (4.39).  r  is  varied  to  represent  gain 
modulation.  For  Eqs.  (4.43)  and  (4.44).  p  is  varied  to  represent  loss  modulation. 
From  their  definitions  in  Eqs.  (4.8)  and  (4.15).  these  parameters  can  be  measured 
experimentally. 

4.4  A  Method  for  Determining  Intrinsic  Cavity  Loss 

The  threshold  for  a  laser  system  is  defined  to  be  the  operating  point  at  which  the 
round  trip  gain  is  equal  to  round  trip  loss.  In  other  words,  the  rate  at  which  the 
signal  grows  inside  the  cavity  is  balanced  by  its  decay  rate  out  of  the  cavity.  In  this 
section  we  derive  an  equation  that  represents  the  threshold  condition  and  describe  an 
experimental  method  to  determine  the  intrinsic  losses  of  a  laser  system. 

Let  Q  represent  the  gain  coefficient  (fractional  gain  per  watt  per  second)  of  the 
laser  and  be  the  pumping  power  at  threshold.  Then  the  threshold  condition  can 
be  written  as 

=  (4.49) 

where  the  cavity  decay  rate  7^.  is  defined  in  Eq.  (4.2).  If  we  let  the  total  cavity  decay 
rate  7-  be  represented  by  the  sum  of  the  decay  rate  due  to  intrinsic  losses  of  the 
laser  and  an  induced  decay  rate  ~iaom  due  to  a  variable  loss  element  such  as  an 
acousto-optic  modulator  (.AOM).  then  we  can  write 

!■:  =  lint  +  ';AO\f  •  (4.  )0) 

The  threshold  condition  then  becomes 

~ltnt  +  lAOM  =  . 
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(4.51) 


If  file  losses  are  small,  we  can  express  Kq.  (I.')!)  in  terms  of  percent  loss.  I  sing 
Ecj.  (4.2)  vve  write 

f,-  =  >,ni  +  ( AO.M  =  th  ■  (4..')2) 

Here  g  represents  the  gain  of  the  laser  in  ptiiTfut  gain  ptr  pass  ptr  watt  of  pumping 
power.  In  terms  of  the  fractional  gain  coefficient  a  and  the  cavity  round  trip  time  T 


we  express  g  as 


g  =  lOOro  . 


(4..i:n 


Since  the  intrinsic  cavity  loss  Enf  is  a  constant  of  the  laser,  we  rewrite  Eq.  (  l.')2)  with 
a.s  a  dependent  variable 


(ao.u  ~  • 


I  {..')4) 


This  is  an  equation  representing  a  straight  line  of  slope  g  and  intercept  Erom 

Eq.  (4.o4).  one  can  determine  the  intrinsic  cavity  loss  by  measuring  the  threshold 
pump  power  Wth  required  at  different  levels  of  induced  loss  (.AO\f-  Indeed  this  was 
the  objective  of  the  second  set  of  experiments  described  in  Chapter  o. 

From  this  same  experimenta  method,  we  can  also  determine  other  parameters  of 
our  laser  system.  Recalling  the  definitions  of  the  laser  pump  parameter  r  in  Eq.  (4.8) 
and  the  normalized  pumping  rate  R  in  Eq.  (4.46),  we  can  write 

^  looijr’ 

where  we  have  also  used  Eqs.  (4.49)  and  (4.53).  Therefore,  once  we  know  g  from 
the  slope  of  Eq.  (4.54).  we  can  determine  the  normalized  pumping  rate  R.  which  is 
another  constant  of  the  system,  for  a  given  pump  power  IV’p.  VVe  apply  this  method 
to  determine  R  for  our  laser  system  in  Chapter  5. 

I  inally.  we  can  appro.ximate  the  stimulated  emission  cross  section  <721  of  the  laser 
gain  medium.  .Near  threshold,  the  gain  coefficient  is  known  as  the  small  signal  gain 
coefficient.  In  terms  of  percent  gain  per  watt  per  pass,  the  small  signal  gain  coefficient 
g  can  be  written 

100/?oV)'72^;3'72i7-2 

g  =  -  -  .  ( 
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I  he  //  torins  represent  efficiencies  of  transferring  energy  from  the  pninp  laser  to  the 
laser  transition.  The  quantum  efficiency  r/o  represents  the  fraction  of  absorherl  [)ho- 
tons  into  the  pump  band  that  actually  participate  in  the  laser  transition.  The  energy 
efficiency  term  rq  represents  the  ratio  of  the  energy  of  a  laser  photon  to  that  of  a 
pump  photon 

h  vi  Ap 

V1-7— =  -r- 

hvp  A/ 

Here  the  subscript  p  refers  to  the  pump  and  /  refers  to  the  laser.  The  fraction  of  the 
actual  pump  beam  that  is  incident  on  the  laser  gain  medium  is  represented  by  ;/j. 
which  accounts  for  reflection  losses,  scattering,  etc.  The  fraction  of  the  incident  light 
actually  absorbed  by  the  gain  medium  is  given  by  The  upper  state  lifetime  r> 
is  defined  in  Eq.  (4.4)  and  .4  is  the  cross  sectional  area  of  the  pump  beam  passing 
through  the  gain  medium.  Finally,  hvp  represents  the  energy  of  a  pump  photon. 
Solving  Eq.  (4. 06)  for  (T21  we  find 


gAhup 

<72i  =  rr~ - ^ - .  (f.^S) 

l00r/ot/lJ?2^/37'2 

.Ml  terms  except  go  in  Eq.  (4.58)  can  be  determined  experimentally,  with  a  little  extra 
work,  from  the  method  for  determining  cavity  losses  outlined  in  this  section.  W’e  will 
use  Eq.  (4.58)  in  Chapter  5  to  approximate  the  stimulated  emission  cross  section  (T21 
for  our  Ti:sapphire  crystal. 

In  this  chapter  we  saw  how  Eqs.  (4.1)  and  (4.3)  can  be  used  to  analyze  the  relax¬ 
ation  oscillations  of  lasers.  We  also  developed  a  method  for  modeling  this  dynamical 
behavior  in  lasers.  We  concluded  the  chapter  with  a  brief  look  at  the  threshold  con¬ 
dition  and  saw  how  the  intrinsic  cavity  loss,  the  normalized  pumping  rate,  and  the 
stimulated  emission  cross  section  can  be  determined  with  experimental  measurements 
at  threshold.  In  the  next  chapter,  w’e  describe  the  experiments  conducted  with  our 
Tirsapphire  la.ser. 
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Chapter  5 


Experimental  Results 


In  the  Ti:sapphire  laser,  population  inversion  evolves  on  a  time  scale  comparable 
to  the  cavity  lifetime.  Such  lasers  are  known  as  class-B  lasers.  Most  solid  state 
lasers  fall  into  this  category.  When  perturbed,  this  class  of  lasers  e.xhibits  a  wide 
variety  of  dynamical  behavior.  These  include  simple  relaxation  oscillations,  spiking 
oscillations,  and  chaotic  behavior.  The  type  of  dynamical  behavior  observed  depends 
on  the  frequency  and  depth  of  perturbations  from  the  steady  state. 

This  chapter  begins  with  a  description  of  the  laser  system  used  to  conduct  the 
experiments.  We  then  describe  the  data  acquisition  procedure,  results,  and  analysis 
of  two  sets  of  experiments.  The  first  set  involved  the  measurement  of  the  relaxation 
oscillations  and  spiking  behavior  exhibited  by  the  Ti:sapphire  laser.  The  second  set 
involves  experiments  that  allow  us  to  determine  the  intrinsic  loss  of  the  cavity  by  an 
independent  method.  From  the  analysis  of  these  experiments,  we  were  able  to  derive 
useful  information  such  as  the  upper  lasing  level  decay  rate,  cavity  losses  and  gain 
per  pass  for  our  laser  system.  .Although  the  results  presented  here  are  specific  to  the 
Ti:sapphire  laser  we  designed  and  constructed,  similar  procedures  and  conclusions 
could  be  applied  to  other  laser  systems  as  well.  Information  from  these  experiments 
would  be  useful  in  other  experiments  on  the  dynamics  of  class-B  lasers. 

5.1  Laser  Description 

The  laser  system  used  for  the  experiments  described  in  this  chapter  was  designed 
according  to  the  methods  developed  in  Chapter  3.  A  diagram  of  the  la.ser  is  shown 
in  Fig.  5.1.  The  laser  cavity  consists  of  two  curved  mirrors,  each  with  radius  of 
curvature  R  =  \0  cm  and  two  flat  mirrors  for  steering  the  long  arm.  .All  four  mirrors 
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Figure  o.l:  Laser  cavity  design  used  for  the  experiments  in  Chapter 

were  coated  for  high  reflectivity  in  the  wavelength  range  6o0  to  !)00  nm  and  hish 
transmittivity  (>8o^^)  for  the  argon  ion  pumping  lines  of  488  and  oi  l.")  mii.  I  he 
coating  was  optimized  for  incident  angles  of  20"  for  p-polarization.  The  mirror  mounts 
were  constructed  of  anodized  aluminum  and  had  three  contact  points  for  adjustment. 
This  mount  design  afforded  precise  translational  and  sensitive  lilt  adjustments  with 
80-pifch  adjusting  screws.  The  mirror  mounts  were  secured  to  a  2-inch  stainless  steel 
optical  breadboard  mounted  on  a  four  pedestal  Boating  optics  table. 

The  2  cm  long  Ti:.sapphire  rod  was  cut  with  Brewster  windows  at  both  ends  and 
had  a  diameter  of  4  mm.  The  rod  was  housed  in  a  water-cooled  copper  jacket  mounted 
on  a  tilt/rotation  stage.  This  stage  allowed  for  rotational  adjustment  about  ail  axes 
of  the  crystal.  The  crystal  was  wrapped  with  thin  flexible  copper  foil  to  ensure  good 
thermal  contact  with  the  cooling  jacket. 

The  dimensions  of  the  cavity  were  determined  by  using  the  stability  requirements 
discussed  in  Section  3.2.  The  size  of  the  laser  w'as  to  be  large  enough  to  accommodate 
sarious  intracavity  elements  and  still  fit  the  confines  of  the  1'  x  2'  breadboard.  The 
actual  length  of  the  long  arm  was  approximately  107  cm  and  the  value  selected  from 
the  stability  diagrams  for  b  was  about  1.5  cm  (see  Figs.  3.2  and  .3.3).  This  gave  a 
total  distance  of  11.5  cm  between  the  two  curved  mirrors  along  the  cavity  axis  and 
4.75  cm  between  the  faces  of  the  crystal  and  the  cvirved  mirrors.  The  physical  round 
trip  length  of  the  cavity  was  118  cm.  The  angle  of  incidence  at  the  curved  mirrors 
was  15.5'.  selected  to  give  maximtim  astigmatic  compensation  in  the  long  arm  of  the 
laser,  as  detailed  in  S<x’tion  3.3. 


liie  pump  boam  was  mode  matched  to  the  fundamental  cavity  mode  with  a 
10  cm  achromatic  lens  positioned  7.9  cm  behind  one  of  the  curved  cavity  mirrors, 
as  described  in  Section  3.4.  The  theoretical  position  of  7.3  cm  from  E(}.  (3.68)  was 
initially  used,  then  the  lens  was  translated  to  optimize  output.  The  mount  for  the  lens 
afforded  two  dimensions  of  translational  freedom  as  well  as  tilt  around  the  vertical 
and  lateral  axes.  The  polarization  of  the  pump  beam  was  rotated  with  a  half  wave 
plate  (polarization  rotator)  anti- reflection  coated  for  the  pump  wavelengths.  The 
vertical  polarization  was  rotated  90'^  into  the  tangential  plane  for  minimal  loss  at  the 
Brewster  window  of  the  crystal. 

A  Brewster's  angle  Faraday  rotator  and  compensating  waveplate  were  placed  in 
the  long  arm  of  the  cavity  to  obtain  unidirectional  lasing.  The  combination  was 
placed  inside  the  cavity  to  permit  clockwise  propagation  of  the  laser  beam  in  the  long 
arm,  as  shown  in  Fig.  5.1.  A  three-plate  birefringent  filter,  inserted  at  Brewster's 
angle  was  used  to  obtain  single  wavelength  performance.  Single  mode  operation  was 
obtained  by  introducing  a  3  mm  thick  etalon  into  the  cavity.  The  etalon  has  about 
30  %  reflectivity  centered  at  780  nm  for  normal  incidence.  The  entire  laser  system 
was  enclosed  in  a  plexiglass  dust  cover  to  improve  stability.  The  pump  beam  was 
steered  through  a  plexiglass  tube  to  minimize  fluctuations  due  to  air  currents  and 
dust  particles  over  the  2  meters  of  propagation  to  the  Ti:sapphire  laser.  The  laser 
operated  stably  in  a  single  longitudinal  mode.  This  was  checked  with  the  help  of  a 
scanning  spectrum  analyzer.  Experiments  were  conducted  with  this  laser  as  described 
in  the  following  sections  of  this  chapter. 


5.2  Experiment  1:  Relaxation  Oscillations 

There  were  two  objectives  in  the  first  set  of  experiments.  The  first  objective  was 
to  determine  the  fluorescence  lifetime,  rj.  of  the  upper  lasing  level  and  the  intrinsic 
cavity  losses,  (percent  loss  per  round  trip),  of  our  Ti:sapphire  laser.  Both  of  these 
cpiantities  were  determined  by  analyzing  the  results  of  the  experiments  in  terms  of 
the  linearized  rate  equation  model  introduced  in  Chapter  4. 

The  second  objective  was  to  see  how  well  the  full  rate  equations  predicted  the 


55 


()l)s»>rve<l  behavior  of  the  laser.  The  theoretical  model  assumed  a  four  level  laser  and 
ignore<l  all  sour<-es  of  noise,  fundamental  as  well  as  deterministic.  It  assumefl  that 
polarization  dynamics  are  fast  compared  to  the  population  and  field  dynamics  and 
that  the  lower  level  decays  almost  instantaneously  to  the  lower  pumping  level. 

The  experimental  method  involved  switching  the  cavity  Q.  or  modulating  the  laser 
cavity  loss.  The  output  waveforms  from  the  laser  were  digitized  and  recorded  as  a 
function  of  time.  .A  series  of  measurements  were  made  of  the  relaxation  oscillations 
at  different  pump  parameters.  From  this  series,  one  particular  waveform  was  selected 
at  random  to  e.xtract  characteristic  laser  parameters.  Once  known,  these  parameters 
were  used  to  predict  the  relaxation  oscillation  behavior  of  the  laser  at  other  pump 
parameters. 


5.2,1  Experimental  Procedure 

In  order  to  perturb  the  laser,  an  acousto-optic  modulator  (.AO.\I)  was  placed  in  the 
long  arm  of  the  laser,  as  shown  in  Fig.  5.2.  The  .AO.M  was  manufactured  by  Intra.Ac- 
tion  Corporation.  It  consists  of  a  quartz  crystal,  anti-reflection  coated  for  minimal 
static  insertion  loss  in  the  700  to  900  nm  wavelength  range,  sandwiched  between 
two  piezoelectric  transducers.  The  transducers  set  up  acoustic  waves  in  the  crystal 
causing  spatially  periodic  changes  in  the  index  of  refraction.  These  periodic  changes 
give  rise  to  a  refractive  index  grating  in  the  path  of  the  beam.  This  grating  causes 
a  certain  fraction  of  the  beam  to  be  diffracted  out  of  the  beam  path  and  acts  as  a 
source  of  intracavity  loss.  The  loss  suffered  by  the  beam  passing  through  the  .AOM 
depends  on  the  angle  of  incidence  of  the  beam  and  the  amplitude  of  the  acoustic  wave. 
■An  SO  .\IHz  acoustic  wave  was  set  up  by  the  driver  (Intra.Action  Corporation  Model 
■MF-SOB).  The  driver  produces  an  SO  MHz  signal  whose  amplitude  is  proportional  to 
the  input  voltage  in  the  range  0  to  1  volt.  The  depth  and  frequency  of  modulation 
(duration  of  acoustic  power  on)  was  controlled  by  pulses  from  a  Krohn-Hite  .Model 
2000  signal  generator  to  the  driver. 

.After  the  .AO.M  was  put  inside  the  cavity,  all  elements  were  adjusted  to  minimize 
the  threshold  of  laser  action.  Laser  threshold  was  found  to  be  at  4.75  watts  from  the 
argon  ion  pump  laser  when  no  signal  was  applied  to  the  .AOM.  The  cavity  was  then 
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Figure  5.2:  Experimental  setup  to  conduct  relaxation  oscillation  experiments  de¬ 
scribed  in  Section  5.2.  The  laser  beam  is  represented  by  the  light  dashed  line  and  the 
pump  beam  by  the  heavy  dashed  line.  The  solid  line  represents  electrical  conn'Ttions. 


Figure  5.5:  .\  timing  diagram  comparing  the  outputs  of  the  signal  generator,  the 
.AOM  driver,  and  the  laser  for  shallow  cavity  loss  modulation.  The  cavity  experienced 
greater  loss  during  the  time  period  that  the  acoustic  wave  was  on. 


loss  modulated.  Both  the  frequency  and  depth  of  loss  modulation  could  be  controlled. 
Cavity  losses  are  related  to  cavity  Q  by 


2ttL 

■ 


(.j.i) 


where  L  is  the  round  trip  length.  A  is  the  lasing  wavelength,  and  /,  is  the  percent 
power  loss  per  round  trip.  It  follows  that  since  the  .AOM  modulates  the  internal  loss, 
it  modulates  the  cavity  Q. 

Output  from  the  laser  w'as  detected  by  a  variable  gain,  high  bandwidth,  fast 
photodetector  using  the  light  output  from  one  of  the  flat  cavity  mirrors.  To  keep  the 
threshold  pump  power  as  low  as  possible,  no  output  coupler  was  used.  This  setup  also 
allowed  us  to  determine  just  the  intrinsic  cavity  loss,  without  having  to  compensate 
of  output  losses.  The  signal  from  the  photodetector  was  monitored  with  a  Tektronics 
246o.A  doO  -MHz  oscilloscope  and  digitized  with  a  LeCroy  TR8837F  transient  recorder. 
The  transient  digitizer  was  housed  in  a  CAM  AC  crate  controlled  by  a  LeCroy  8901. A 
GPIB  crate  controller.  Both  the  transient  digitizer  and  the  .AO.M  driver  were  triggered 
from  the  same  source.  .At  the  end  of  each  measurement,  the  recorded  waveforms  were 
transferred  to  an  IB.VI  PC/.AT  for  further  analysis.  A  full  schematic  of  the  e.xperiment 
is  shown  in  Fig.  5.2.  .A  timing  diagram  comparing  the  output  of  the  signal  generator, 
the  .AO.M  driver,  and  the  laser  as  a  function  of  time  is  shown  in  Fig.  5.3.  Data  were 
taken  and  evaluated  for  two  independent  experiments  to  determine  the  parameters 
that  govern  the  rate  equations  of  Section  4.1. 


5.2.2  Results 

For  the  transient  dynamics  experiments,  the  cavity  Q  was  modulated  at  two  different 
depths.  Deep  modulation  resulted  in  the  laser  being  completely  off  during  the  pe- 
riocl  of  modulation.  This  means  that  the  cavity  loss  and  .AO.M  induced  loss  exceeded 
round  trip  gain.  Shallow  modulation  resulted  in  the  laser  operating  continuously,  but 
jumping  between  different  "cavities."  one  with  a  high  Q  and  the  other  with  a  lower 
Q.  For  the  sake  of  distinguishing  the  response  of  the  laser  to  the  loss  modulation  in 
these  two  cases,  we  refer  to  the  transient  oscillations  in  the  case  of  deep  modulation 


as  .^p/knnj  ox  illations  and  those  in  the  case  of  shallow  modulation  as  irlaxatiori  oscil- 
I  he  output  waveforms  were  recorded  for  different  pump  parameters.  At  each 
operating  point,  we  recorded  waveforms  from  a  single  sweej)  and  from  an  average  of 
ten  sweeps. 

One  particular  waveform  averaged  over  ten  sweeps,  chosen  at  random,  was  ana¬ 
lyzed  to  determine  all  the  characteristic  parameters  describing  the  relaxation  oscilla¬ 
tion  equations  of  Section  4.2.  The  procedure  for  determining  T2  and  the  corresponding 
uncertainty  is  as  follows.  From  Eq.  (4.27).  v\-e  see  that  the  envelop  of  transient  os¬ 
cillation  is  a  pure  exponential  with  a  decay  rate  3.  This  exponential  decay  rate  is 
related  to  the  upper  state  decay  rate  */2  via  Eq.  (4.30).  Therefore,  knowing  the  pump 
parameter  r  and  measuring  3  allows  us  to  determine  '2.  The  upper  state  lifetime  is 
related  to  the  upper  state  decay  rate  by  the  expression 


To  meeisure  3.  we  let  fi{t)  represent  the  peak  amplitudes  of  the  oscillations  from 
the  steady  state.  These  peaks  form  an  exponentially  decaying  curve  according  to  the 
expression 

pit)  =  .  (0.3) 


Here  represents  the  amplitude  of  the  initial  perturbation  and  3  is  the  decay  rate. 
The  positions  of  the  peaks  of  the  oscillation  versus  elapsed  time  were  fit  to  an  expo¬ 
nential  decay  curve  to  determine  3  in  Eq.  (5.3).  This  measured  value  of  3  was  then 
used  along  with  the  pump  parameter  of  the  waveform  being  analyzed  to  determine  the 
upper  state  decay  rate  according  to  Eq.  (4.30).  We  then  determine  T2  using  Eq.  (5.2). 

The  uncertainty  in  T2  was  estimated  from  the  uncertainty  in  i.  Since  the  peaks 
of  the  recorder!  waveform  correspond  to  a  specific  number  of  counts.  .V.  we  assume 
a  Poissonian  distribution  for  the  counts  and  take  the  statistical  uncertainty  in  the 
peaks  to  be 


.V 


v/T 


=  /v. 


(5.4) 


The  steady  state  value  of  the  waveform  vvas  determined  by  taking  the  time  average 
of  the  las  le  half  of  the  waveform,  where  the  laser  had  rewhed  steadv  state.  I'he 


uncertainty  in  the  steady  state  value  was  also  estimated.  The  statistical  uncertainty 
in  time  separation  between  successive  peaks  was  simply  the  root  mean  squared  of 
the  tleviations  of  the  time  period  between  peaks.  If  we  take  the  natural  logarithm 
of  Ec|.  (5.3).  we  can  e.xpress  3  in  terms  of  the  amplitudes  of  the  peaks  and  the  time 
periods,  T, 


3  = 


In  —  In 
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Using  the  laws  of  error  propagation  [36],  we  can  approximate  the  uncertainty.  A3,  in 
3  to  be 


A3 
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Here  A/i  represents  the  averaged  uncertainty  of  the  peaks  of  the  waveform.  Afi.,  rep¬ 
resents  the  uncertainty  in  the  steady  state  value  of  the  waveform,  and  AT  represents 
the  uncertainty  in  the  time  separation  between  peaks. 

Using  r  =  1.39  for  the  pump  parameter  of  the  waveform  being  analyzed  and  the 
measured  value  of  3  in 


2J 


12  = 


(0.7) 


we  find  the  value  *^or  the  upper  state  decay  rate  to  be 


72  =  288  ±  10  KHz  . 


(5.8) 


The  upper  state  lifetime  T2  is  then  given  by 


72  =  3.47  ±  .13^isec  . 


(•5.9) 


To  find  U.  the  intrinsic  cavity  loss,  we  first  determine  the  relaxation  oscillation 
frequency  u,'r.  from  the  meastired  waveform  using  the  average  time  period  betwefui 
successive  peaks  Tn,, 

27- 

(5.10) 
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From  the  measured  value  of  we  are  able  to  determine  the  cavity  loss  rate.  *■...  by 
using  Eq.  (4.16)  into  (4.29)  to  obtain 

,2  .2  .  .12 
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(5.11) 
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I  he  uncertainty  in  the  oscillation  frequency  was  determined  from  the  uncertainty  in 
the  time  periods  according  to 

A^-.  Ar 
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The  uncertainty  in  ye  is  determined  from  the  relation 
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Using  the  measured  values  for  3.  "2  and  r  in  Eq.  (5.11),  '/e  is  found  to  be 


(5.l:l) 


(5.14) 


ye  =  9.18  ±  .05  MHz  . 


(5.15) 


Using  Eq.  (5.15)  into  Eq.  (4.2)  we  obtain  a  cavity  loss  of 

4  =  (.3.6  ±  .2)  %  per  pass 


(5.1b) 


for  our  1 18  cm  cavity.  This  figure  includes  the  static  insertion  loss  of  the  .-XOM  which, 
according  to  the  specifications,  carries  a  loss  of  about  0.25  %. 

Once  the  values  for  the  upper  state  decay  rate  y2  and  the  cavity  decay  rate  y,-  are 
determined,  a  rela.xation  oscillation  frequency  can  be  predicted  for  all  pump  parame¬ 
ters  and  com.pared  with  the  rela.xation  oscillation  frequency  measured  experimentally 
from  the  recorded  waveforms.  Figure  5.4  shows  the  results  of  such  a  comparison 
for  deep  modulation.  All  experimental  data  for  vVr  were  extracted  from  single  sweep 
records.  Figure  5.6  shows  the  same  comparison  except  that  the  experimental  data 
wen’  extracted  from  waveforms  averaged  over  ten  sweeps.  Error  bars  for  the  single 
sweep  results  represent  the  root  mean  squared  deviations  in  the  period  of  oscillation. 
1  he  error  bars  for  the  ten  sweep  results  represent  the  largest  variation  between  the 
single  sweep  measurement  and  the  ten  sweep  measurement  for  a  given  pump  parame¬ 
ter.  For  completeness,  a  plot  of  the  output  power  of  the  laser,  measured  by  the  same 
detector  used  for  detecting  the  waveform  signal,  as  a  function  of  pump  parameter 
is  shown  in  Fig.  5.5.  Figures  5. 7-5. 9  show  the  same  comparison  as  Figs.  5.4 -5.6  for 
shallow  modulation  of  the  cavity. 
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Figure  o.  1:  Variation  of  oscillation  frequency  o-v  with  pump  parameter  r.  Fxpenmen- 
tal  points  are  from  single  sweep  waveforms  recorded  with  deep  cavity  loss  modulation. 
The  continuous  line  represents  the  theoretical  curve  from  Eq.  (4.29).  Parameters 
and  *.2  were  extracted  from  a  single  waveform  averaged  over  ten  sweeps  at  the  pump 
parameter  r  =  1.39. 


PUMP  PARAMETER  r 

Figure  o.o:  Variation  of  laser  output  power  versus  pump  parameter  r  for  deep  mod 
Illation.  The  line  represents  theoretical  output  characteristics. 
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PUMP  PARAMETER  r 

Figure  5.(i;  Variation  of  oscillation  frequency  a-v  with  pump  parameter  r.  l-.xp<  i  iin-  ii 
tal  points  are  from  waveforms  averaged  over  ten  sweeps.  The  data  point  with  the 
cross  indicates  the  waveform  used  to  extract  parameters  for  the  theoretical  curves  in 
Figs.  0.4.  0.6.  0.7  and  5.9. 


PUMP  PARAMETER  r 

Figure  Variation  of  oscillation  frequency  with  pump  parameter  r  lor  slialiow 
cavity  loss  modulation.  Experimental  points  are  from  single  sweep  waveforms.  The 
continuous  curve  represents  theoretical  predictions  based  on  Eq.  (4.29).  Parameters 
and  *2  were  extracted  from  a  single  deep  modiilation  waveform  averaged  over  ten 
sweeps  for  pump  parameter  r  =  1.39. 
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PUMP  PARAMETER  r 

Figure  0.8;  Variation  of  laser  output  power  as  a  function  of  pump  para  met  tr 
shallow  modulation.  The  line  represents  theoretical  output  characteristics. 


PUMP  PARAMETER  r 

Figure  ■).!):  Variation  of  oscillation  frequency  with  pump  parameter  r  for  .slial 
low  cavity  modulation.  E.xperimental  points  are  from  waveforms  averaged  over 
ten  sweeps.  The  continuous  curve  represents  the  theoretical  predictions  based  on 
Eq.  (4.29).  Parameters  and  *2  were  extracted  from  a  single  deep  modtilalion 
waveform  averaged  over  ten  sweeps  for  pump  parameter  r  =  1.39. 
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1  hi-^  ('xperiiiuMit  was  repeated  to  explore  the  role  of  noise,  intrinsic  or  extrinsic,  in 
the  experiment.  I  he  same  experimental  procedure  was  used,  with  the  exception  that 
only  single  sweep  waveforms  were  recorded.  .\11  parameters  were  then  extracted  from 
a  single  waveform  as  opposed  to  a  waveform  averagefl  over  ten  sweeps.  I  he  result  for 
the  upper  state  lifetime  was 

r2  =  3.48  ±  .lS//sec  .  (o.lT) 

and  the  cavity  loss  was  determined  to  be 

f,-  =  (3.6  ±.4)  7c  per  pass  .  I 

The  plots  in  Figs.  5.10  through  5.13  show  the  observed  and  predicted  spiking  and 
relaxation  oscillation  frecpiency  and  the  corresponding  output  power  as  a  function  of 
pump  parameter.  No  error  bars  are  shown  in  Fig.  5.12  since  in  many  cases  there  were 
only  two  peaks  on  the  waveform.  This  single  data  poin’^  did  not  allow  an  estimate  of 
error  in  the  period  or  frequency  of  oscillation. 

5.2.3  Discussion  of  Results 

The  close  agreement  between  the  two  data  collection  runs  for  the  values  for  the  up¬ 
per  state  lifetime  T2.  and  the  intrinsic  cavity  loss  suggests  that  noise  did  not 
significantly  effect  the  experimental  measurement  of  these  quantities,  e.xcept  per¬ 
haps  at  lower  operating  points.  The  value  for  the  upper  state  lifetime  can  be  com¬ 
pared  with  other  measurements  whose  primary  objective  was  to  measure  this  value 
precisely.  Spectroscopic  and  lasing  characteristics  of  Ti:sapphire  have  been  of  con¬ 
siderable  interest  ever  since  its  potential  as  a  broadly  tunable  solid  state  laser  was 
recognized  [1.  37.  38].  Experiments  to  determine  the  fluorescence  lifetime  of  the  Ti 
ioi.  in  .AI2O3  have  been  carried  out  by  monitoring  the  fluorescence  decay  after  exciting 
a  crystal  with  a  pulse  of  radiation  from  a  dye  laser  at  various  wavelengths  [1.  39.  40]. 
The  results  of  various  experiments  are  consistent  and  show  a  temperature-dependent 
upper  state  lifetime.  The  lifetime  has  a  maximum  value  of  approximately  3.85  psec 
below  200  K.  It  then  decreases  rapidly  in  the  temperature  range  of  200  to  400  K  to 
a  value  less  than  1  //sec. 
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PUMP  PARAMETER  r 

l  igure  5.10:  Variation  of  oscillation  frequency  with  pump  parameter  /  .  Kxpfi; 
mental  points  are  from  single  sweep  waveforms  recorded  with  deep  cavity  loss  mod¬ 
ulation.  The  continuous  line  represents  the  theoretical  curve  from  Eq.  (4.29).  The 
data  point  with  the  cross  indicates  the  waveform  used  to  extract  parameters  and 
‘,2  for  Figs.  5.10  and  5.12. 


Figure  5.1  1:  Variation  of  laser  output  power  as  a  function  of  pump  parameter  r  for 
deep  modulation  in  the  second  experiment.  The  line  represents  theoretical  output 
<  haracteristics. 


66 


PUMP  PARAMETER  r 


Figure  5.12:  Variation  of  oscillation  frequency  av  with  pump  parameter  r  for  shallow 
cavity  loss  modulation.  E.xperimental  points  are  from  single  sweep  waveforms.  The 
continuous  line  represents  theoretical  predictions  based  on  Eq.  (4.29).  Parameters 
and  *2  were  extracted  from  a  single  sweep  waveform  for  pump  parameter  r  =  1 .5S. 


PUMP  PARAMETER  r 

Figure  5.14:  \’ariation  of  la.ser  output  power  as  a  function  of  pump  parameter  r  for 
'hallow  modulation  in  the  second  experiment.  The  line  represents  theoretical  output 
characteristics. 


In  our  experiments,  the  crystal  was  cooled  by  circulating  water  in  a  jacket  sur¬ 
rounding  the  crystal.  Good  thermal  contact  between  the  jacket  and  the  crystal  wa.s 
achieved  by  wrapping  the  crystal  with  thin,  flexible,  copper  foil.  This  cooling  kept 
the  crystal  at  a  temperature  of  approximately  280  K.  Extracting  a  value  for  the  upper 
state  lifetime  from  the  published  results  [1. 39.  40]  for  this  temperature  gives  a  lifetime 
of  3.20  to  3.3o  /<sec.  Our  results  are  consistent  with  this  value,  within  experimental 
error.  This  agreement  gives  support  to  our  experimental  method  for  determining  the 
upper  state  lifetime  T2  by  analyzing  the  relaxation  oscillations.  Further,  this  agree¬ 
ment  suggests  that  the  assumptions  made  in  deriving  the  equations  governing  the 
relaxation  oscillations  in  Section  4.2  were  indeed  valid. 

The  value  of  the  intrinsic  cavity  loss  G  is  unique  to  our  laser  system.  Therefore  it 
can  only  be  compared  with  a  measurement  of  the  cavity  loss  by  a  method  independent 
of  the  relaxation  oscillations.  This  measurement  was  the  primary  objective  of  the 
.second  set  of  experiments,  the  results  of  which  are  described  in  Section  5.3. 

To  see  how  well  the  rate  equation  model  is  able  to  predict  the  observed  behavior  of 
the  laser,  we  compare  the  experimental  results  with  the  theoretical  predictions  based 
on  Eq.  (4.29)  in  Figs.  5.4.  5.6.  5.7,  5.9.  5.10.  and  5.12.  Such  a  comparison  shows  that 
the  general  behavior  of  the  relaxation  oscillation  frequency  follows  the  theoretical 
prediction  at  low  pump  parameters  (r  <  1.5).  At  higher  pump  parameters  however, 
noticeable  deviations  from  the  predicted  behavior  appear.  This  same  behavior  is 
observed  in  the  output  power  of  the  laser,  as  shown  in  Figs.  5.5.  5.8.  5.11.  and  5.13. 
The  output  power  should  increase  linearly  with  increased  pump  power,  as  depicted 
by  the  continuous  line  in  these  figures.  This  deviation  of  laser  output  from  a  linear 
dependence  on  the  pump  parameter  has  been  observed  in  other  experiments  [22, 
41],  This  departure  from  the  expected  output  power  is  probably  due  to  temperature 
changes  inside  the  mode  volume  of  the  crystal  at  higher  pump  powers.  .\s  previously 
noted,  the  upper  state  lifetime  is  temperature  dependent,  with  a  steep  slope  in 
the  room  temperature  operating  regime.  This  means  that  a  small  change  in  the 
temperature  inside  the  crystal  can  result  in  a  significant  change  in  the  upper  state 
lifetime.  Recalling  the  small  gain  coefficient  given  in  Eq.  (4.56).  we  see  that  the  gain 
of  the  laser  has  a  linear  dependence  on  the  upper  state  lifetime.  .A  shorter  upper 


state  lifetime  reduces  the  gain  per  pass  inside  the  laser  and  hence  would  decrease  the 
output.  This  would  have  the  same  effect  as  lowering  the  pump  parameter  r.  since  the 
threshold  for  lower  gain  requires  greater  pumping. 

The  effect  that  a  change  in  the  upper  state  lifetime  would  have  on  the  relaxation 
oscillation  frequency  can  be  determined  by  analyzing  Eq.  (4.29).  The  dominant  fre¬ 
quency  is  ul'o,  found  in  the  first  term  of  the  square  root  and  defined  in  Eq.  (4.16). 
Since  r  decreases  with  a  decreasing  upper  state  lifetime,  we  would  anticmate  that  the 
rela.xation  oscillation  frequency  would  decrease  as  the  temperature  insi  *he  crystal 
goes  up.  This  is  what  is  seen  in  Figs.  5.4.  5.6,  5.7,  5.9.  5.10.  and  5.12.  .\t  lower  pump 
parameters,  in  the  range  of  r  =  1.0  to  1.5,  the  observed  relaxation  oscillation  fre¬ 
quency  agrees  with  the  predictions  of  Eq.  (4.29)  derived  in  Section  4.2.  This  suggests 
that  the  simple  laser  rate  equations  (4.1)  and  (4.3)  describe  the  dynamics  of  the  laser 
system  well  for  lower  pump  parameters,  where  thermal  effects  can  be  ignored. 

5.2.4  Results  of  Numerical  Modeling 

We  also  analyzed  the  data  obtained  from  our  experiments  in  terms  of  a  model  that 
describes  the  dynamics  of  the  system.  This  comparison  between  theory  and  experi¬ 
ment  is  more  detailed  than  that  based  solely  on  Eq.  (4.29)  and  provides  a  much  more 
stringent  check  of  the  vahdity  of  the  simple  rate  equation  model.  In  Section  4.3.  a 
method  for  modeling  the  laser  rate  equations  (4.1)  and  (4.3)  of  Section  4.1  was  pre¬ 
sented.  For  a  loss  modulated  system,  the  appropriate  equations  of  motion  were  given 
by  E'.qs.  (4.43)  and  (4.44).  These  equations  were  solved  numerically  with  experimen¬ 
tally  measured  parameters  using  the  Runge-Kutta  method  of  algorithms.  Computer- 
generated  waveforms  were  then  produced  based  on  these  solutions  for  comparison 
with  the  recorded  waveforms.  Using  this  technique,  we  modeled  several  waveforms 
selected  at  random.  The  method  in  which  the  computerized  model  "laser"  was  used 
to  generated  the  various  waveforms  is  as  follows. 

For  deep  modulation,  the  model  ef|uations  were  solved  starting  at  /  =  0  with  the 
number  of  atoms  in  the  upper  lasing  state  at  the  scaled  stationary  state  value  of  p 
as  defined  in  Eq.  (4.48).  Using  the  results  from  our  first  experiment  for  *  i  (5.8)  and 
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(5.19) 


*  .  (5.15)  into  K(i.  (  t.  l5).  the  equation  for  the  modulating  parameter,  we  find 

/>  =  —  a:  .31.9  . 

“2 

The  normalized  pumping  rate  R  \va.s  set  as  defined  by  Eq.  (4.46)  using  the  pump 
parameter  of  ♦he  waveform  being  modeled  and  the  value  of  p  given  in  Eq.  (5.19).  W’e 
started  with  the  laser  "off"  initially.  For  this,  the  value  for  the  number  of  photons 
in  the  cavity  was  first  estimated  to  be  10“'*^.  based  on  the  value  of  the  saturation 
photon  number  q,  defined  in  Eq.  (4.33).  The  actual  number  of  photons  in  the  cavity 
mode  is  indeterminant  because  of  spontaneous  emission  noise.  The  initial  number  of 
photons  was  adjusted  to  get  the  model  laser  to  initiate  buildup  at  approximately  the 
same  time  as  the  recorded  waveform.  With  the  parameters  initialized  as  described, 
the  waveform  generated  by  the  model  represents  the  turn  on  of  a  laser  recovering 
from  deep  modulation. 

Comparisons  of  the  measured  and  model-generated  waveforms  for  deep  cavity 
loss  modulation  are  shown  in  Figs.  o.  14-5. 1 7.  The  model-generated  waveforms  are 
represented  by  the  continuous  curve  and  the  measured  waveform  by  the  dotted  curve. 
Figure  5.14  compares  a  single  sweep  waveform  at  pump  parameter  r  =  1.29  with  the 
model  predictions  and  Fig.  5.15  compares  a  waveform  averaged  over  ten  sweeps  with 
the  model  predictions  at  the  same  pump  parameter.  Figure  5.16  compares  a  single 
sweep  waveform  at  pump  parameter  r  =  1.39  with  the  model  predictions  and  Fig.  5.17 
compares  a  waveform  averaged  over  ten  sweeps  at  the  same  pump  parameter. 

For  shallow  modulation,  the  initial  values  for  the  model  laser  were  chosen  in 
the  same  manner  as  for  deep  modulation,  except  that  the  cavity  photon  number 
was  set  to  its  stationary  state  value  (r  —  1).  as  defined  by  Eq.  (4.47).  The  loss 
parameter  p  was  then  modulated  at  a  frequency  corresponding  to  the  modulating 
frequency  of  the  .AOM  in  the  actual  experiment.  The  model  was  allowed  to  "run"  long 
enough  to  eliminate  transient  effects  of  the  initial  conditions.  Since  the  exact  depth 
of  cavity  loss  modulation  was  not  known  for  the  relaxation  oscillation  experiments, 
the  depth  of  modulation  for  p  was  determined  by  varying  the  fractional  change  in  p 
until  the  steady  state  values  of  the  recorded  and  modeled  waveforms  overlapped.  The 
actual  depth  of  modulation  for  the  cavity  loss  was  then  estimated  from  this  fractional 
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Figure  5.14:  A  comparison  of  recorded  and  modeled  waveforms  representing  the  spik¬ 
ing  oscillations  of  the  Ti:sapphire  laser  with  deep  loss  modulation.  The  solid  curve 
depicts  the  modeled  waveform  generated  by  the  model  described  in  Section  4.3.  The 
dotted  curve  represents  the  single  sweep  waveform  recorded  at  r  =  1.29, 


Figure  5.15:  A  comparison  of  recorded  and  modeled  waveforms  representing  the 
spiking  oscillations  from  deep  loss  modulation.  The  solid  curve  depicts  the  model¬ 
generated  waveform  and  the  dotted  curve  represents  the  recorded  waveform  averaged 
over  ten  sweeps  at  r  =  1.29. 
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f  igure  *1.1^:  A  comparison  of  recorded  and  modeled  waveforms  representing  tii*'  -pik¬ 
ing  oscillations  from  deep  loss  modulation  at  pump  parameter  r  =  The  solid 

curve  depicts  the  waveform  generated  by  the  model  and  the  dotted  curve  represents 
the  recorded  single  sweep  waveform. 


1  iitiirf  ").1T:  A  comparison  of  recorded  and  modeled  waveforms  representing  tin' 
spiking  oscillations  from  deep  loss  modulation.  The  solid  curve  depicts  the  model- 
generated  waveform  and  the  dotted  curve  represents  the  recorded  waveform  a\eraged 
over  ten  sweeps  at  r  =  1.39. 
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rhauHt'  in  p.  With  the  parameters  for  the  model  deterniined  as  described  here,  the 
getierated  \va\ef()rnis  represent  a  laser  system  being  distnrberl  wifli  shallow  cavity 
loss  modulation. 

Figures  o.lS-o.'il  show  the  comparison  between  the  model  predictions  and  the 
measured  waveforms  for  shallow  cavity  loss  modulation.  The  model  generated  wave¬ 
forms  were  adjusted  on  time  scale  to  compensate  for  the  trigger  delay  in  the  n*corded 
waveforms.  The  depth  of  loss  modulation  for  all  modeled  waveforms  was  about  0.3  'T. 
Figure  5.18  shows  a  comparison  of  a  single  sweep  waveform  recorded  at  pump  param¬ 
eter  r  =  1.28  with  the  model  prediction.  waveform  averaged  ovet  ten  sweeps  at 
the  same  pump  parameter  compared  with  the  model  prediction  is  shown  in  h'ig.  5.1!). 
Figures  5.20  and  5.21  compare  the  waveforms  recorded  at  another  pump  parameter 
of  value  r  =  1.50  for  a  single  sweep  and  averaged  over  ten  sweeps  with  the  respective 
model  predictions. 

By  comparing  the  recorded  waveforms  with  the  model-generated  waveforms  we 
conclude  that  the  model  describes  the  qualitative  behavior  of  the  laser  from  loss 
modulation  well.  Quantitative  comparisons  between  the  model  and  the  actual  laser 
behavior  are  limited  by  a  number  of  factors.  First  is  the  fact  that  the  recording  of 
the  actual  laser  behavior  by  the  transient  digitizer  is  limited  by  the  speed  of  the 
digitizer.  .-\  comparison  between  the  signal  rise  time  on  the  oscilloscope  and  the 
rise  time  displayed  on  the  recorded  waveform  indicated  tha'  the  bandwidth  of  the 
transient  recorder  was  less  than  that  of  the  detector.  Second,  since  the  model  uses 
experimentally  measured  parameters,  uncertainty  in  the  measured  values  will  result 
in  uncertainty  in  the  modeled  waveforms. 

We  can  also  derive  conclusions  based  on  the  differences  between  the  single  swr-ep 
wa\eforms  and  those  averaged  over  ten  sweeps.  For  the  deep  modulation  waveforms, 
the  averaged  waveforms  consistently  have  shallower  initial  peaks  or  "spikes"  than 
the  waveforms  recorded  from  single  sweeps.  For  shallow  modulation  th.-  difference 
is  not  as  pronounced.  Slight  fluctuations  in  the  time  to  reach  this  Hrst  peak  would 
result  in  an  overall  decrease  in  the  amplitude  of  the  spike  as  it  is  averaged.  This 
suggests  that  the  rise  time  to  the  first  peak  fluctuates  more  from  a  cavity  with  def'p 
loss  modulation  than  for  a  cavity  with  shallow  loss  modulation.  1  h('  source  of  this 


Figure  o.l8:  A  comparison  of  recorded  and  modeled  waveforms  representing  the  re¬ 
laxation  oscillations  of  the  Ti:sapphire  laser  with  shallow  loss  modulation.  The  solid 
curve  depicts  the  modeled  waveform  generated  by  the  model  described  in  Section  1.3. 
The  dotted  curve  represents  the  single  sweep  waveform  recorded  at  r  =  1.28.  Depth 
uf  loss  modulation  for  the  model  was  0.3  %. 


TIME  (Msec) 

l  igure  ■).l!);  .\  comparison  of  recorded  and  modeled  waveforms  representing  the  re¬ 
laxation  oscillations  from  shallow  loss  modulation.  The  solid  curve  depicts  the  model¬ 
generated  waveform  and  the  dotted  curve  represents  the  recorded  waveform  averaged 
over  ten  sweeps  at  r  =  1.28. 
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Figure  5.20:  A  comparison  of  recorded  and  modeled  waveforms  representing  the  re¬ 
laxation  oscillations  from  shallow  loss  modulation  at  pump  parameter  r  =  1.50.  The 
solid  curve  depicts  the  waveform  generated  by  the  model  and  the  dotted  curve  rep¬ 
resents  the  recorded  single  sweep  waveform.  Depth  of  loss  modulation  wa.'  D.  l  /'. 


Figure  5,21;  comparison  of  recorded  and  modeled  waveforms  repre.senting  the  re¬ 
laxation  cjscillations  from  shallow  loss  modulation.  The  solid  curve  depicts  the  model- 
generated  waveform.  The  dotted  curve  represents  the  recordcxl  waveform  averaged 
over  ten  sweeps  at  r  =  1.50. 


slioht  tt'in[)oral  Huctuation  for  tiie  <l«'p  modulated  cavity  is  (juautum  noise  since  tiu' 
laser  i"  otf  tlnring  the  time  the  cavity  suffers  large  loss.  The  buildup  time  to  tin* 
first  [)eak  is  dependent  upon  the  number  of  photons  inside  the  fundamental  mode  of 
the  cavity.  I  he  buildup  of  the  laser  radiation  starts  once  a  population  inversion  has 
developed  and  a  spontaneously  emitted  photon  pas.ses  through  the  gain  medium.  I  he 
average  time  for  laser  buildup  decreases  as  the  number  of  photons  present  initially  in 
the  cavity  increases.  The  study  of  this  quantum  noise  and  its  deterministic  effect  on 
the  rise  time  of  a  laser  signal  offers  an  excellent  opportunity  for  further  research. 

5.3  Experiment  2:  Cavity  Losses 

In  Section  o.’2.  we  described  an  experiment  to  measure  the  intrinsic  loss  of  a  laser 
system  by  studying  its  rela.xation  oscillations.  In  this  section,  we  describe  aiu.)ther 
experiment  that  was  carried  out  to  determine  the  intrinsic  cavity  losses,  fhis  rnc'thod 
also  allows  us  to  determine  the  normalized  pumping  rate  R  of  our  laser  as  defined  by 
Eq.  (  l.lfi)  and  estimate  the  stimulated  emission  cross  section  rr^i  of  the  Tiisapphire 
crystal.  The  theoretical  background  for  this  experiment  is  given  in  Section  J.  1.  I’o 
carry  out  this  experiment,  the  loss  offered  by  the  acousto-optic  modulator  used  in  the 
experiment  was  calibrated.  The  experimental  procedure  consisted  of  measuring  the 
threshold  pump  power  as  a  function  of  the  .AOM  loss. 

5.3.1  Experimental  Setup:  AOM  Calibration 

Our  objective  in  calibrating  the  .\0.\I  loss  was  to  determine  the  percent  loss  out  of 
t  h<'  primary  Iteam  as  a  function  of  input  voltage  to  the  driver  of  the  .\0M.  The  .AOM 
was  placed  at  normal  incidence  immediately  outside  the  I  i:sap[)hire  laser  cavity,  flu' 
laser  was  then  oi)erated  at  a  certain  pump  power  above  threshold,  fhe  light  output 
from  one  of  the  cavity  mirrors  was  allowed  to  pass  through  the  .\0.\f.  f  lie  zero  ord('r 
beam  in  the  transmitted  light  illuminated  a  photodeterMor.  fhe  detector  was  [)laced 
far  enough  away  from  the  .AO.M  to  ensure  sufficient  separation  Ix'tween  th<'  zero  aiul 
higher  order  diffracted  lieatns  so  that  they  would  not  be  detecterl.  Ihe  signal  trom 
the  photodetector  was  displayed  on  an  o.scilloscope  and  digitized  and  recorded  b\  the 
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I  raii'K'ut  (liiiilizt’i'.  I'lu'  source  of  DC  voltage  to  the  AOM  ilriver  was  a  Pow('r  Designs. 
liK  .  Miiilel  tiOoOA  DC  power  supply.  Flie  carrit'r  levr'l  (hiasi  oii  tin’  driver  was  sei  to 
z<‘ro.  1  he  DC  \oltag<'  to  the  driver  was  iiureineiited  iti  steps  id  l(j(J  iii\  in  the  ranuf' 
Ironi  0  to  1  \  .  \\ a\('torins  wc're  reconh'd  at  eacfi  level  of  DC  input  to  the  dii\er. 
D)  luitiiniize  the  (effect  of  fluctuations  in  the  pump  and  the  rCsapphire  laser,  tin' 
waveforrtis  were  averaged  over  lUl)  sweeps.  I  he  recorded  waveforms  were  then  further 
time  averaged  over  the  2  m.s  duration  of  the  swet'p  to  determine  an  average  value  for 
the  light  intensity.  .\ny  contril)ution  from  the  detector  offset  was  subtracted  from 
the  measured  mean  signal  producing  a  number  where  .r  represents  the  a|)[)lied 
voltage  to  the  driver  and  .V  represe'tits  the  average  intensity  of  the  transmitted  beam. 
Percent  loss  for  the  ap|)lied  voltage  was  determined  according  to  tlie  relation 

k)ss  induced  by  .\0M  =  — — -  x  100  .  (  ').20) 

•^0 

The  data  points  showing  9?  loss  as  a  function  of  DC  volts  applied  to  the  driver  ar(' 
shown  in  Fig.  ■■).22.  T  he  straight  line  is  the  best  Ht  to  the  data  points.  The  slope 
of  this  line  gives  the  relationship  between  loss  induced  out  of  the  primary  beam  and 
voltage  applied  to  the  .AO.M  driver.  This  relationship  was  determined  to  be 

loss  =  0.0l2l9f/mV  .  (^.dl  ) 

I  he  uncertainty  in  this  value  is  less  than  O.lo'/f  of  the  measured  value. 

5.3.2  Experimental  Setup:  Cavity  Losses 

Onie  the  .\0.\I  was  calibrated,  it  was  placed  at  normal  incidence  in  the  long  arm  ol 
the  (  avity.  .\11  intracavity  ehunents  were  adjusted  to  minimize'  the'  thn'shold.  which 
was  determined  to  be  at  l.To  watts  of  pump  power.  Output  from  the  last>r  was 
monitored  by  a  high  gain  detector  from  reflections  off  one  of  the  Drewster  windows 
e)f  the'  crvstal  as  slujwn  in  Fig.  1  his  light  was  filtered  with  a  colored  glas'^ 

filter  and  focused  onto  the  elete'ctor.  I  he  signal  from  tht'  df'tector  was  monitore'd 
with  an  oscilloscope.  The  threshold  e'oiulition  was  ilcfined  as  the  point  at  whieh  the' 
fluctuations  in  the  signal  were  equal  to  the  average  of  the  signal.  1  he  intt'iisity  of  tlie 


Figure  5.22:  Foss  calibration  curve  for  the  IntraAction  Corp.  acousto-optic  nuutuiauu 
used  in  the  relaxation  oscillation  and  the  cavity  loss  cxp«*rinienls.  Data  were  taken 
at  normal  incidence  using  output  trom  the  Fi'.sapphire  laser. 

pttrnp  power  was  measured  using  a  photodetector  monitoring  th(‘  reHections  ol  i)um[) 
light  from  the  cavity  curved  mirror,  as  shown  in  Fig.  5.2.F  The  rellectf’d  li^iit  was 
attenuated  through  a  series  of  neutral  density  filters. 

Voltage  applied  to  tfie  .\0.\I  driver  was  incremented  in  steps  of  approximately 
;50  m\'.  The  range  of  applied  voltage  was  from  zero  to  a  maximum  of  approximately 
200  m\  .  .\t  the  level  of  loss  corresponding  to  200  mV  to  the  driver,  the  pum[)  power 
was  insufficient  to  overcome  the  .VOM  induced  loss.  .\t  each  increment  of  voltage  to 
the  .VOM  driver,  the  pump  power  was  increa.sed  until  the  laser  threshold  condition  was 
met.  For  each  value  of  voltage  applied  (mV)  to  the  .VO.M.  threshold  pump  power  was 
re(  (;rfled  frcjin  the  output  of  the  detector  in  volts.  I  his  [)roc<'ss  was  repeated  thri'c 
tiiiK's  in  succession.  Kach  measurement  was  convertefl  into  data  jioints  as  follows, 
file  \()ltage  appli<'<l  to  the  VO.Vl  driver  was  converted  into  i)ercent  loss  basivl  cui  the 
formula  given  in  Fip  (5.21).  1  he  signal  from  the  photodetector  was  converte<l  into 

watt>  of  pump  power  by  using  th.e  formula 

Watts  =  [().!)  X  signal  voltage]  x  2.0  (  >.22' 


I  bis  formula  ua.s  determined  Iw  calil)rating  theontjuit  voltagi'  from  the*  photodetcctor 
against  the  out|)Ut  [lower  readings  from  the  meter  on  the  power  sipijily  of  the  pump 


Lix  r.  lilt'  \i)lrai>(‘  from  the  <l(neetor  was  reconied  at  two  values  of  pump  power  (  s:  d 
au<l  ID  watts)  as  indicated  on  tlu'  meter  on  the  power  supply.  Ihis  voltai>e  was 
then  reducf'd  hy  t  h(“  \'oltage  reading  ohtaiiied  from  the  di'tector  when  the  iiu  idont 
Ix'am  was  blocked  with  a  Ix'am  dmnp  and  ordy  background  light  was  dett'cted.  rhi> 
reduction  represented  appro.ximatelv  \i)V(  of  the  ■■unl)lo<ke(r‘  valut'.  Tht'  resulting 
voltage  readings  corresponding  to  the  two  pump  power  s<'ttings  wer<'  comi)ared  with 
the  power  meter  reatliugs  to  establisli  a  scale  factor  of  ^  ‘2. 


5.3.3  Results 


-Ml  the  data  points  of  ‘X  loss  with  threshold  pump  ]>ower  were  iIkmi  plottt'd  and  litted 
with  a  line  of  minimum  regression.  The  results  are  shown  in  Fig.  o.LM.  I  he  data 
points  are  shown  for  only  two  of  the  three  runs  to  minintize  overcrowding  of  the  data 
points.  The  straight  line  depicts  the  best  fit  to  all  data  points.  Tlx*  ititercejjt  of 
the  straight  line  on  the  ( ,\o\f  axis  gives  the  value  of  the  intrinsic  cavity  loss,  in 
accordance  with  Eci.  (4.54).  The  measured  value  for  the  intrinsic  cavity  loss  i  from 
these  measurements  is 


(,  =  (2.95  ±  .35)'X, 


(5.23) 


In  order  to  determine  the  nortnalized  pumping  rate  R  deiined  by  Ftp  (  l.lb).  we 
measured  the  slope  of  the  line  of  regression.  .\s  described  in  Section  4.4.  the  slopi' 
represents  <y.  the  percent  gain  per  pass  through  the  cavity  |)er  watt  of  pump  |)ower. 
Ihe  nieasured  value  of  g  was 


g  =:  (0.7DT  ±  0.01  1 )  'X  [ter  watt  per  pass 


( 5.2  1 1 


I  "iing  this  value  in  Fq.  (1.55)  with  the  upper  state  decay  rate  dett-rmined  from  the 
first  experiment  we  find 

R  =  [7.30  i  .38]  X  U  ,,  .  (  5.25) 

1  his  experimc'ut  also  allows  us  to  estimate  the  gain  cross  section  tjf  the-  me<lium 
with  the  help  of  Fq.  (  1.50).  I  lu*  eHi<  ieticies  in  Fq.  (4.50)  were  »'stimated  as  follows. 
Ihe  rpiantum  efficiency  could  not  be  measured  directly,  but  we  <'stimat('d  it  to  be 


^  0.8 


I  5.2fi ) 
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Figure  u.23:  1  he  experimental  setup  for  conducting  the  cavity  loss  experiment  using 
the  AOM  to  \ary  cavity  loss.  The  dashed  lines  represent  optical  paths  and  the  -(did 
lines  represent  electrical  paths. 


THRESHOLD  PUMP  POWER  (wotts) 


1  luurt'  \'ariali(uis  of  threshold  pum()  power  with  intracavity  .\().\1  indui  c-d  !o>> 

1  he  data  points  are  shown  for  two  t)f  the  three  runs.  The  straight  line  represents  the 
best  lit  to  all  data  [)oints. 


ii.i'ci!  wii  |)it'\  ii>u>  t'Xjx'ririK'nts  '1.  10.  12].  1  lie  (MU'rgv  etticiency  l|^  was  (leteniiiiicil 

in>in  !  l.’iT  I  usiii^  jOO  nm  as  the  average  wav('|(Migth  i)f  the  leirnp  ptuaoiis  and 
!>()  nil!  a.'  the  ax'erage  wavelength  of  a  laser  [ihotiin 

t/i  O.fi  f  .  I  -i.d  I  I 

1  he  primary  loss  in  the  transmission  of  the  pump  bt'am  to  the  crystal  was  from 
reflections  off  the  curved  cavity  mirror  prior  to  the  crystal.  Measurements  of  thi' 
reflected  beam  compared  with  the  output  from  puntp  la.ser  showed  a  consistent  loss 
of  28'/.  This  gives  an  efficiency  t/o 


t/2  =  0.72  .  !  ’).2>i 

Single  pass  absorption  measurements  of  the  crystal  showed  that  apitroximately  >>.:?'/  of 
the  pump  beam  was  absorbed  by  the  crystal  over  a  pumping  range  of  O.o  to  7.0  watts. 
I  his  gives 

V3  =  O.Sd  .  i.').20) 

lo  determine  .1.  the  average  cross  sectional  area  of  the  pump  beam  inside  the  crystal, 
we  refer  to  the  mode  matching  problem  described  in  Section  2. 1  and  the  characteristics 
of  the  spot  size  at  the  beam  waist  given  in  Section  2.2.  For  our  cavity,  the  spot  size 
at  the  beam  waist  was  estimated  to  be  2.2  microns  (see  Figs.  2.1  an<l  2.-')).  Since  the 
crystal  is  2  cm  long,  we  take  the  average  spot  size  to  be  at  a  position  -'lO  mm  inside  tlu' 
crystal  Ir  =  ."tO  mm).  To  determine  the  spot  size  at  this  location,  we  use  K(|s.  (2.1  1) 
and  (  2.40 )  to  obtain 

tr.,,,  =  10  microns  .  (')..20) 

iiiving  an  average  cross  sectional  area  of  the  pump  beam 

.4  ^  •'»  X  10“’  cm"’  .  I  Tdl  ) 

1  sing  t  hese  valties  in  Fq.  (  l..')8).  wo  estimate  t  he  t  he  st  imulated  emission  cross  si'ct  ion 
to  be 

fT^i  1.1  X  10“''*  crn’  .  ('t.22i 


SI 


5.3.4  Discussion  of  Results 


lilt'  dost'  Hitreeiueiit  in  the  measured  values  for  the  intrinsic  cavity  loss  I.,  bt'lwt't'ii 
this  t'X[)eriuH'nt  and  the  relaxation  oscillation  experiments  suggests  that  Uoth  t'xpt'r- 
ituental  techniqties  are  reliable.  Indeed,  a  [u'ocedure  ha.seil  on  variable  loss  similar  to 
that  used  in  the  second  set  of  experiments  has  been  used  by  Sanchez  r/  al.  to  measure 
the  intrinsic  loss  of  a  Ti;sapphire  la.ser  [41].  However,  the  technique  of  determining 
cavity  lo.s.ses  by  analyzing  the  relaxation  oscillations  of  a  laser  is  a  new  technique  that 
can  be  easily  ap{)lied  to  other  solid  stale  lasers. 

for  our  results,  the  slightly  higher  value  for  the  intrinsic  cavity  loss  from  the  Hrst 
experiment  is  probably  due  to  an  increase  in  thermal  losses  with  increased  pumping 
power.  The  values  of  loss  from  the  first  set  of  experiments  were  measured  frotn  data 
obtained  while  the  crystal  was  being  pumped  with  approximately  6.5  watts  from  the 
pump  laser  (r  =  1.49  and  IT,;,  =  4.75  watts).  Comparing  this  pump  power  with  the 
pump  power  used  in  the  second  set  of  experiments  (see  Fig.  5.24),  we  see  that  the 
data  from  the  second  experiment  was  for  the  most  part  obtained  using  lower  pump 
powers. 

Schulz  gives  an  excellent  account  of  the  thermal  effects  in  the  i'itsapphire  crystal 
in  [22].  He  estimates  that  between  30  and  45  percent  of  the  power  absorbed  by  the 
crystal  is  deposited  as  heat,  resulting  in  thermal  lensing  and  a  change  in  the  imlex  of 
refraction  inside  the  crystal.  Both  of  these  effects  have  the  result  of  ■■detuning"  the 
cavity  or  increasing  loss. 

I  he  gain  cross  section  has  been  measured  using  several  different  experimental 
t<'chniques  fl.  38.  40.  11.  12.  13].  The  measured  values  range  from  6.5  x  10“-’*'  end 
to  1.5  X  If)-'*'*  cm^  depending  on  technique  and  wavelength  used  in  the  measurement. 
Our  result  agrees  well  with  these  mea.su rements. 

In  this  chapter  we  have  described  in  detail  the  laser  system  used  to  conduct  two 
different  experiments.  Fhe  simple  laser  rate  equation  approximations  developed  in 
Chapter  1  seem  to  describe  well  the  dynamical  behavior  of  the  Ti:sapphire  relaxation 
oscillations.  The  cavity  loss  experiments  confirmed  that  th('  intritisic  cavity  losses  (  an 
be  determined  from  relaxation  o.scillations.  Both  sets  of  experiments  open  ui>  many 
opportunities  for  studying  the  dynamics  of  ri:saf)phire  lasers.  Sotiu'  of  th«'s(’  will  !)e 


<ii'(  u<-c(l  ill  ( 'ha|)tt-'r  i . 
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Chapter  6 


Polarization  Properties  of  Gaussian  Beams 


Our  treatment  of’  laser  beams  in  Chapter  2  was  within  the  framework  of  the  scalar 
wave  equation  (2.7)  and  the  paraxial  approximations  (2.2.3)  and  (2.24).  For  most 
applications  which  do  not  involve  the  polarization  properties  of  laser  beams,  this 
framework  is  quite  adequate.  Indeed,  the  cavity  design  and  experiments  described 
in  this  thesis  were  developed  within  this  framework.  When  describing  the  polariza¬ 
tion  properties  of  laser  beams  however,  a  scalar  representation  is  inadequate.  Kven 
for  a  linearly  polarized  laser  beam,  a  scalar  description  of  a  finite  cross  section  laser 
beam  is  inconsistent  with  .VlaxwelFs  equations  (2.1)-(2.4).  The  transverse  nature 
of  the  electromagnetic  field  expressed  by  Eqs.  (2.1)  and  (2.3)  implies  that  the  spa¬ 
tial  variation  of  the  field  in  directions  transverse  to  the  direction  of  propagation  is 
coupled  to  the  polarization  properties  of  the  field.  Thus  it  is  well  known  that  spa¬ 
tial  variation  of  the  field  in  the  direction  of  polarization  gives  rise  to  a  longitudinal 
field  component  [44].  This  coupling  of  the  transverse  spatial  variation  of  nonplanar 
wavefronts  to  polarization  was  investigated  in  an  interesting  paper  by  Fainman  and 
Shamir  [l-'i].  They  analyzed  the  cross  polarization  in  <i  spherical  wavefront  from  a 
point  source.  They  also  recorded  experimentally  the  cross  polarization  of  a  linearly 
polarized  fundamental  Gaussian  beam  passing  through  a  pin  hole.  Simon.  Sudar- 
■ihan  and  Mukunda  [44l  used  an  analysis  of  Maxwell's  equations  in  the  relativistic 
frcuit  form  to  establish  the  polarization  properties  of  the  fundamental  Gaussian  laser 
beams.  Several  other  approaches  to  this  problem  have  been  discus.sed.  including  po¬ 
tential  formulation  of  Maxwell's  equations  [46.  47]  and  a  power  series  expansion  of 
electromagnetic  fields  [48].  In  this  chapter,  we  present  another  approach  to  this  prob¬ 
lem.  Our  approach  is  simpler  and  more  direct  and  allows  us  to  csstablish  the  general 
polarization  structure  of  paraxial  Hermite^Gaussian  modes  of  a  laser.  We  also  pres('nt 


al  exideiicp  for  this  struct ure  tisin|>,  the  linearly  polarizt'd  beam  from  an 
argon  ion  laser. 

6.1  Paraxial  Solutions  to  Maxwell’s  Equations 

In  Chapter  2.  we  derived  the  paraxial  wave  ecpiat ion  (2.26)  and  developed  a  solution 
representing  the  basi<'  characterist ics  of'  Gaussian  beams.  In  this  section,  we  deter¬ 
mine  the  form  of  paraxial  solutions  to  Maxwell's  etjuations  which  allow  us  to  belter 
represent  the  polarization  properties  of  laser  beams.  For  (piasi-  monochromatic  fields 
propagating  in  the  c-direction.  we  write  the  electric  and  magnetic  fields  of  the  wave 
as 

E{rJ)=  E(r)c‘'*’'-— ■'>  =  (/TdrlG  +  •  '''-H 

B(r.O=  =  («,(r)G  +  .  (6.2) 

Here  E(r)  and  B(r)  describe  the  transverse  spatial  profile  of  the  beam  and  f_..  and 
(  ire  unit  vectors  along  the  .r.  y,  and  r-axes  respectively.  For  the  fields  itt  t(|s.  |6.1 ) 
and  (6.2)  Maxwell's  equations  in  free  space  become 


/A-/-:3(r)-F  V-E(r) 

=  0  . 

(6.2) 

ikfr^  X  E(r)  -F  V  X  E(r) 

==  /A-cBtr)  . 

( 6. 1 1 

ikB-sir)  +  V  •  B(r) 

=  0  . 

( b..-!) 

ik'ts  X  B(r)  +  V  X  B(r) 

=  -t-E(r)  . 

(6.6) 

c 


Frotn  these  equations  we  see  that  each  cartesian  component  of  the  electric  and  mag- 
netif  field  satisfies  the  scalar  wave  equation  (2.7).  hor  paraxial  beam  like  solutions, 
the  inequalities  defining  the  paraxial  approximation  Kqs.  (2.22)  and  (2.2!)  hold  tot- 
each  field  component  and  each  component  satishes  the  paraxial  wave  eqiiat  ioti  (2.26). 
These  componenls  are  not  arbitrary.  They  are  cotipled  via  Fqs.  (6.2)  (6.6).  In  the 
paraxial  approximation  of  Eqs.  (2.22)  and  (2.2t),  F<js.  (fi.-l)  and  (b.d)  allow  u.s  to 
express  the  longit tidinal  fiehl  components  in  terms  of  the  transverse  field  components 
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where  we  have  used  the  inequality  from  Kq.  (2. If:!'  in  arriving  at  these  eciuatious. 

I  sing  Ec|.(6. 1)  we  can  express  the  magnetic  Held  components  in  terms  of  the 
electric  field  components  as 
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Similarly,  using  Eq.  (6.7)  and  the  paraxial  wave  equation  (2.26)  in  Eq.  (6.6).  we  can 
express  the  electric  field  components  in  terms  of  the  magnetic  field  components 
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In  writing  Eqs.  (6.7)  and  (6.8),  we  have  kept  terms  up  to  \/{kw]^,  where  ir  is  some 
characteristic  length  scale  associated  with  the  transverse  beam  profile.  .\n  inspection 
of  these  eciuations  shows  that  the  electric  and  magnetic  field  components  can  be 
expressed  in  terms  of  two  solutions  of  the  paraxial  wave  equation  (2.26).  Let  us 
denote  these  solutions  of  the  paraxial  scalar  wave  equation  by  /(r)  and  ,v(r).  then 
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One  can  easily  check  that  Eqs.  (6.1o)-(h.20)  satisfy  the  paraxial  Maxwell's  equa¬ 
tions  (6.3)-(6.6)  up  to  terms  of  order  l/(/‘W’)^.  We  now  consider  some  explicit  exam¬ 
ples  of  these  solutions. 


6.2  Linearly  Polarized  Electromagnetic  Waves 

From  Eqs.  (6.io)-(6.20)  we  see  that  finite  cross  section  electromagnetic  waves,  in 
general,  have  a  longitudinal  field  component.  In  this  sense  it  is  not  possible  to  have 
pure  transverse  electromagnetic  beams.  However,  since  the  longitudinal  component 
is  smaller  by  a  factor  of  1/ktc  compared  to  the  transverse  fields,  it  is  possible  to  have 
beams  that  have  dominant  transverse  polarization.  We  now  examine  the  form  of  a 
linearly  polarized  electromagnetic  beam. 

Without  loss  of  generality,  we  take  the  direction  of  dominant  polarization  to  be 
the  j-direction.  Then,  choosing 
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we  find  Eqs.  (6.15)-(6.17) 

take  the  form 
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where  c',n-i(r)  is  given  by  Eq.  (2. -58). 

In  writing  Eq.  (6.22)  we  have  dropped  terms  involving  second  order  derivatives 
of  c’mn(r)  since  compared  to  the  first  order  terms  they  are  smaller  by  the  factor 
l/(A-u’)^.  The  leading  term  in  the  cross-polarization  term  is  smaller  by  the 

factor  l/(fcti-)^  compared  to  the  £’|-component.  It  is  kept  here  because  if  the  beam 


l)y  Kfi".  ((i.'J'i)  (6.21)  is  [)asse(l  through  a  crossfvl  polarizer,  tlu'ii  is  tlio 

only  torrn  that  is  t rarisrnitte<l.  I  he  longitudinal  term  is  smaller  by  the  lactor  [/htr 
compared  to  the  dominant  transverse  component. 

The  magnetic  Held  components  correspoiuling  to  Eqs.  (H.2‘2)-(t).2  1)  are  obtained 
by  the  relations 
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L  sing  the  properties  of  Hermite  polynomials  we  can  write  down  explicit  expressions 
for  E.2""‘*(r)  and  as 
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These  ecpiations  explicitly  indicate  the  relative  magnitudes  of  various  fi(“ld  compo¬ 
nents.  A  more  quantitative  measure  of  the  relative  strengths  of  various  field  compo¬ 
nents  is  obtained  by  comparing  the  powers  associated  with  different  field  components. 
Let  the  total  power  associated  with  the  beam  be  P^.  Then  the  power  associated  with 
different  electric  field  components  is  found  to  be 
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It  follows  from  these  equations  that  the  power  associated  with  different  field  compo¬ 
nents  is  a  constant  of  propagation,  that  is.  it  is  independent  of  c.  Equations  (b.dl ) 
(6.dd)  together  with  Eqs.  (b.28)-(6.30)  show  that  an  electromagnetic  beam  predoiTi- 
inantly  polarized  in  the  .r-direction  has  a  small  cross-polarization  <ompon<’nt  in  llu' 
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(/-ilirecriuii  in  addition  to  a  small  longitudinal  component.  1  he  longitudinal  field 
component  is  smaller  hy  a  factor  of  l/Auc  and  the  cross-polarization  Held  component 
is  smaller  hv  a  factor  of  compared  to  the  dominant  polarization  comjjonent. 

From  these  considerations  it  is  clear  that  a  pure  transverse  linearly  polarizer!  electro¬ 
magnetic  beam  is  the  geometrical  optics  limit  (Aue  —*  oc)  of  Kqs.  (6.281  and  (6.2!)). 
In  general,  for  finite  cross-section  beams,  both  the  cros.s-polarization  and  lotigitiidinal 
field  components  must  be  kept  for  consistency  with  .Ma.xweU's  erpiations. 

The  presence  of  longitudinal  field  components  is  also  recjuired  for  a  correct  de¬ 
scription  of  energy  flow  in  the  beam.  Energy  flow  in  a  beam  is  described  in  terms  of 
rays,  which  are  curves  along  which  energy  is  transported.  The  tangent  to  these  curves 
at  a  given  point  indicates  the  direction  of  energy  flow  of  the  Poynting  vector.  For  an 
electromagnetic  beam,  rays  should  converge  as  they  approach  a  focal  region  (beam 
waist)  and  diverge  as  they  leave  the  focal  region.  The  Poynting  vector  averaged  over 
an  optical  cycle  is 

S  =  ;ye<,c^3fe(E  X  B)  .  (6. 21) 

Using  the  fields  given  by  Eqs.  (6.28)-(6.30)  and  certain  recursion  relations  for  Hermite 
polynomials,  we  find  that  the  Poynting  vector  for  the  beam  can  be  written  as 

S  =  -eoC^|TTOnr^|Umn(*’)|^  ‘^f-l  +  ‘^^2  +  ■  (fi.-Io) 

The  quantity  before  the  square  brackets  is  simply  the  beam  intensity  (watts/m‘^). 
The  vector  inside  the  square  brackets  denotes  the  direction  of  energy  flow.  If  we 
recall  that  the  radius  of  curvature  is  negative  for  a  converging  wave  and  positive  tor 
a  diverging  wave,  it  follows  from  Eq.  (6.25)  that  for  a  beam  approaching  the  beam 
waist  (focal  region),  energy  flow  occurs  tow^ard  the  axis  and  away  from  the  ;  axis  as 
it  leaves  the  focal  region.  If  the  longitudinal  component  of  the  fields  is  ignored,  then 
the  energy  flow  occurs  only  along  the  c-axis  and  such  a  beam  is  not  only  inconsistent 
with  Maxwell's  equations,  but  also  does  not  provide  a  correct  description  of  focusing 
properties  of  laser  beams.  From  E(j.  (6.2-5)  we  find  the  equation  for  the  family  of  rays 
is 

P  =  P'\J{~ho)^  +  I  (6.26) 
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wlicrc  p  =  y.r-  4-  //-  and  p,,  is  the  distance  of  the  ray  from  the  f)eam  axis  at  tlie  beam 
waist  ;  =  0. 

I  h('  cross-polarization  component  is  even  smaller  than  the  ionsitiulinal  c(jmpo- 
nent.  If  polarization  properties  of  the  wave  are  not  of  interest,  it  may  he  ignored. 
For  a  correct  description  of  polarization  properties  however,  the  cross- polarization 
component  must  be  kept.  For  example,  if  the  beam  passes  through  a  lint'ar  polar 
izer  whose  axis  is  cros.sed  with  respect  to  the  dominant  direction  of  polarization,  the 
cross-polarization  component  is  the  dominant  component  in  the  transmitted  beam. 

It  is  al.so  interesting  to  compare  the  distribution  of  fields  in  a  plane  transverse  to 
the  direction  of  propagation.  In  general,  the  transverse  distribution  of  fields  evolves 
during  propagation.  It  is  interesting  to  note  that  this  evolution  does  not  involve 
a  change  in  the  energy  associated  with  the  field.  This  follows  from  Fqs.  (b.dl) 
(6.33)  where  we  showed  that  the  power  associated  with  each  field  component  remains 
constant  during  propagation.  Expressions  for  the  transverse  distribution  of  fields  are 
complicated  in  general.  In  the  far  zone  however,  the  fields  take  the  form 
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These  field  components  lead  to  the  following  distributions  for  the  intensities 
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(6.37) 

(6.3S) 

(6.3b) 


(6.10) 

(6.41) 

(6.42) 


higures  6.1.  6.2.  and  6.3  show  these  intensity  distributions  for  TF.Moo  mode  (laussian 
beams.  For  each  of  these  figures,  the  c-axis  scale  is  arbitrary  and  not  related  betwevn 
figures.  In  the  next  section  we  look  at  the  form  of  a  circularly  polarized  wave. 


Figure  6.2:  .l-dimensional  profile  of  the  intensity  distribution  of  a  (■r()s.>  pulariA-d 

TEMoo  mode  Gaussian  beam  as  defined  by  Eq.  (6.41).  The  c-axis  scale  is  arbitrary. 


Figure  6.:i:  A  ;i-(limensional  profile  of  the  intensity  distribution  of  the  longitufiinal 
component  of  a  TEM(x)  mode  Gaussian  beam  as  defined  by  Ecp  (6.42).  I  he  r-axis 
scale  is  arbitrary. 


6.3  Circularly  Polarized  Electromagnetic  Waves 


For  a  circularly  polarized  wave,  we  take 

fir)  =  ti’mn(r)  and 

Then  various  field  components  are 
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In  the  far  zone,  these  equations  lead  to  the  following  transverse  intensity  distributions 
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6.4  Experimental  Observations 


Eixperimental  observations  of  the  intensity  distributions  for  a  linearly  polarized  laser 
beam  were  made  using  an  argon  ion  laser  operating  at  488  nm  with  a  TEMto  mode. 


Our  (>l)ie(  ri\e  was  to  coiiipare  the  linearly  polarized  intensity  prolile  ot  the  laser  beam 
aiiven  Ity  Kq.  (b.lO)  with  its  cross  polarized  intensity  profile  given  in  K(|.  (b.ll  I. 

1  he  intensity  profile  of  t  h<'  laser  beam  was  photographed  using  a  Cohu  Inc.  l''l() 
series  monochrome  solid-state  CCl)  camera.  The  picture  was  recorded  using  an  IHM 
PC  and  Beamcode  6.1  software.  .\  picture  of  the  linearly  polarized  beaiti  was  tak(Mi 
with  the  laser  operating  at  less  than  10  milliwatts  (minimum  scale  on  the  digital 
power  meter)  and  the  beam  attenuated  through  a  neutral  density  filter  with  optical 
density  (OD)  2.  Pictures  of  the  beam  intensity  profile  are  reproduced  in  Figs.  6.1 
and  6.0.  Figure  6.4  shows  a  d-dimensional  contour  of  the  beam  intensity  and  Fig.  6..") 
shows  a  2-dimensional  planview  of  the  beam  intensity.  To  observe  the  cross  polarizcvl 
intensity  profile,  the  laser  was  operated  at  approximately  1 10  milliwatts  and  the  beam 
was  transmitted  through  two  dichroic  sheet  polarizers.  The  polarizers  were  oriented 
so  that  their  polarization  axes  were  orthogonal  to  each  other.  Pictures  of  the'  cross 
polarized  beam  intensity  profile  are  reproduced  in  Figs.  6.6  and  6.7.  Figure  6.6  shows 
the  d-dimensional  pattern  and  Fdg.  6.7  shows  the  2-dimensiona)  planview.  I'he  linear 
and  cross  polarized  intensity  distributions  of  Eqs.  (6.40)  and  (6.41)  can  be  clearly 
seen  from  these  pictures.  Due  to  the  fact  that  the  dichroic  materials  polarize  using 
absorption,  the  relative  intensity  magnitudes  of  the  beams  photographed  in  Figs.  6.4- 
6.7  are  undetermined. 

In  this  chapter  we  have  discussed  the  polarization  properties  of  (Jaussian  beams 
by  developing  paraxial  solutions  to  Maxwell's  equations.  We  saw  that  the  field  com¬ 
ponents  of  these  solutions  were  not  independent,  but  coupled  by  Eqs.  (6.d)-(6.6). 
1  his  results  in  the  presence  of  cross  polarized  and  longitudinal  fields  in  the  transverse 
direc  tion.  The  presence  of  the  cross  polarized  field  was  observed  experimentally. 


Figure  6.4:  A  4-dimensioiial  intensity  profile  of  a  linearly  polarized  i  EM(X)  laser  iM-ani. 
The  picture  was  taken  from  an  argon  ion  laser  operating  at  488  nm. 


Figure  6.5;  A  2-dimensional  plain  iew  of  the  intensity  profile  of  a  linearly  polarized 
TEMoo  laser  beam.  Ihe  picture  was  taken  from  an  argon  ion  laser  operating  at 
488  nm. 


Figure  6. 6:  A  :}-dimensional  iiileiisity  pro(ile  of  a  cross  polarized  IE  Moo  last*!’  l>eain. 
The  picture  was  taken  from  an  argon  ion  laser  operating  at  488  nrn. 


Figure  6.7:  A  "i-dimensional  planview  of  the  intensity  profile  of  a  cross  polarized 
TEMoo  laser  beam.  The  picture  was  taken  from  an  argon  ion  laser  operating  at 
488  nm. 


Chapter  7 


Conclusions 


In  this  thesis  we  have  reviewed  the  characteristics  of  laser  propagation  and  used 
those  characteristics  to  develop  a  method  for  designing  a  laser  system.  The  unique 
characteristics  of  the  Ti:sapphire  laser  make  it  an  excellent  choice  for  a  variety  of 
research  applications.  The  techniques  described  in  this  thesis  can  be  used  to  de¬ 
sign  a  Ti;sapphire  or  other  type  of  laser  system  for  specific  research  applications. 
These  techniques  address  the  requirements  of  stability,  astigmatism  compensation, 
and  mode  matching  of  the  pump  and  cavity  modes. 

VVe  also  studied  some  features  of  the  dynamical  behavior  of  Ti:sapphire  lasers.  10 
our  knowledge,  this  is  the  first  study  of  the  relaxation  oscillations  of  a  liisapphire 
laser.  Tsing  rate  equation  theory,  we  developed  a  method  for  measuring  the  upper 
state  lifetime  and  intrinsic  cavity  loss  analyzing  these  oscillations  in  the  Ti:sapphire 
la.ser.  Our  measured  values  for  these  parameters  are  consistent  and  in  agretnnent 
with  the  values  of  the  upper  state  lifetime  and  the  intrinsic  cavity  loss  of  our  laser 
obtained  by  independent  methods.  Our  method  for  determining  these  parameters 
from  the  study  of  relaxation  oscillations  presents  a  new  and  simple  technique  which 
can  be  applied  to  other  laser  systems  as  well. 

Ifv  also  using  rate  equation  theory,  we  developed  a  method  for  modeling  the 
dynamical  behavior  of  lasers.  This  method  scales  the  laser  rate  etjuations  in  such  a 
way  as  to  allow  numerical  solutions  from  experimentally  measured  parameters. 

\Ve  described  an  experiment  for  measuring  the  intrinsic  cavity  loss  by  inducing 
variable  cavity  loss  and  measuring  the  pump  power  required  to  reach  laser  threshold. 
This  experimental  procedure  allows  us  to  also  estimate  the  normalized  pumping  rate 
R.  the  small  signal  gain  coefficient  g.  and  the  stimulated  emission  cross  section  cr^i 
of  anv  laser  svstem. 


1  Inall}'  we  discusserl  the  interesting  Held  of  beani-like  vector  solutions  to  Maxwell'^ 
ecjuat  ions.  These  solutions  characterize  the  polarization  propert  ies  of  (laussiau  beam-. 
We  alsi)  saw  experimental  e\  i deuce  of  t  he  cross  p(»larizat  ion  intensity  proHh'  predi<  tt'd 
by  our  solutions. 

Tlie  field  of  laser  dynamics  is  very  rich  and  the  Ti:sapphire  laser  j^rovides  an  excel¬ 
lent  system  for  studying  many  aspects  of  dynamical  behavior.  By  vaiying  the  d(’[)th 
and  frequency  of  cavity  gain  or  loss  modulation,  an  investigation  into  chaotic  dynam¬ 
ics.  deep  spiking  oscillations,  and  generation  of  specific  waveforms  can  be  conducted. 
Studies  of  directic>nal  mode  competition  and  bistability  can  also  be  carried  out  on 
a  laser  system  similar  to  the  one  (lescribed  in  this  thesis.  Other  interesting  o[)[)or- 
tunities  for  research  include  first  passage  time  studies  in  the  scaling  regime  [id.  oDi. 
antiphase  studies,  and  other  multi-mode  phenomena. 
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